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We present the critical theory of a number of zero temperature phase transitions of quantum 
antiferromagnets and interacting boson systems in two dimensions. The most important example 
is the transition of the 5* = 1/2 square lattice antiferromagnet between the Neel state (which 
breaks spin rotation invariance) and the paramagnetic valence bond solid (which preserves spin 
rotation invariance but breaks lattice symmetries). We show that these two states are separated by 
a second-order quantum phase transition. The critical theory is not expressed in terms of the order 
parameters characterizing either state (as would be the case in Landau-Ginzburg- Wilson theory), 
but involves fractionalized degrees of freedom and an emergent, topological, global conservation 
law. A closely related theory describes the superfluid-insulator transition of bosons at half-filling 
on a square lattice, in which the insulator has a bond density wave order. Similar considerations 
are shown to apply to transitions of antiferromagnets between the valence bond solid and the Z2 
spin liquid: the critical theory has deconfined excitations interacting with an emergent U(l) gauge 
force. We comment on the broader implications of our results for the study of quantum criticality 
in correlated electron systems. 
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point. If the transition is second order, there is inter- 
esting universal singular behavior that is manifested in 
many physical quantities. According to the prevalent 
paradigm largely due to Landau and Ginzburg^, these 
universal critical singularities are associated with long 
wavelength low energy fluctuations of the order param- 
eter degree of freedom. When combined with general 
renormalization group ideas^, this notion provides the 
sophisticated Landau- Ginzburg- Wilson (LGW) theoreti- 
cal framework for thinking about critical phenomena in 
various diverse contexts. Specifically, static critical prop- 
erties at non-zero temperature are supposed to be deter- 
mined from effective models in which all modes other 
than the order parameter have been eliminated. Simi- 
larly, for dynamical critical properties, the only degrees of 
freedom that purportedly need be retained are the order 
parameter and at most a few additional "hydrodynamic" 
modes having slow relaxation times due to conservation 
laws. 

Recent years have seen much interest in the study of 
zero temperature phase transition phenomena in corre- 
lated many body systems. Unlike their thermal counter- 
part, such transitions are often driven by quantum fluc- 
tuation effects and are hence known as 'quantum phase 
transitions''^. Indeed, it has been proposed that proxim- 
ity to quantum critical points (QCPs) separating two dis- 
tinct phases is responsible for the anomalous properties 
of some interesting correlated materials such as, for in- 
stance, the cuprate superconductors. Theoretically, the 
LGW paradigm has thus far provided the basic frame- 
work to examine quantum critical phenomena as well. 
In particular, the critical modes specific to a quantum 
critical point are presumed to be the long distance, long 
time fluctuations of the order parameter, described in a 
continuum field theory. 

In the last few years some interesting and tantaliz- 
ing evidence has emerged that points toward the failure 
of the LGW paradigm at certain quantum phase tran- 
sitions. First, there are numerical calculations'*'^ that 
see a direct second order quantum phase transition be- 
tween two phases with different broken symmetry char- 
acterized by two apparently independent order param- 
eters. A LGW description of the competition between 
such two kinds of orders would then generically predict 
either a first-order transition, or an intermediate region 
of coexistence where both orders simultaneously exist, or 
an intermediate region with neither order. A direct sec- 
ond order transition between these two broken symmetry 
phases would seem to require fine-tuning to a 'multicrit- 
icaP point. Are the numerics managing to achieve this 
'fine-tuning' or is the LGW paradigm simply invalid? 

Second, there have been a number of fascinating ex- 
periments probing the onset of magnetic long range or- 
der in a class of rare-earth inter-metallics known as the 
heavy fermion metals^ii. Remarkably, the behavior right 
at the quantum transition between the magnetic and non- 
magnetic metallic phases is usually very different from 
that of a Fermi liquid. Furthermore, such behavior is in 



severe disagreement with expectations based on LGW 
analyses. Specifically, theories associating the critical 
singularities with fluctuations of the natural magnetic or- 
der parameter in a metallic environment seem to have a 
hard time explaining the observed non- Fermi liquid phe- 
nomena. Once again it appears that more than the ob- 
vious possibly happens at some quantum critical points. 

In this paper we demonstrate and study various spe- 
cific examples of quantum phase transitions which vio- 
late the LGW paradigm. We will show that in a number 
of different quantum transitions, the natural field the- 
oretic description of the critical singularities is not in 
terms of the order parameter field(s) that describe the 
bulk phases but in terms of some new 'emergent' degrees 
of freedom that are specific to the critical point. These 
new degrees of freedom may be thought of 'fractional' 
quantum number particles that interact with each other 
through an emergent gauge force in a sense made pre- 
cise below. Laughlin has previously argued for fraction- 
alization at quantum critical points on phenomenological 
grounds^. A non-technical overview of our results has ap- 
peared previously.'' 

We note, in passing, that there are already numer- 
ous well-documented examples of the breakdown of the 
LGW paradigm in quantum systems in one dimension^ 
However, these rely rather crucially on the description of 
various states in terms of the harmonic phase degrees of 
freedom of the Tomonaga-Luttinger liquid, and do not 
have any direct generalization to higher dimensions. 

For the most part in this paper, we will study phase 
transitions in two dimensional quantum magnetism. 
These may also be fruitfully viewed from a different point 
of view as representing transitions of interacting bosons 
on a lattice at commensurate density. Quantum magnets 
provide a particularly useful laboratory to develop and 
test ideas on the theory of quantum phase transitions. 
Consider a quantum system of spin S = \/2 moments Sr 
on a two dimensional square lattice (r = {x,y)) with the 
Hamiltonian 

H ^ J^Sr- Sr' + ... (1.1) 

(rr') 

The ellipses represent other short ranged interactions 
that may be tuned to drive various zero temperature 
phase transitions. We assume J > 0, i.e antiferromag- 
netic interactions. Later we will consider various gener- 
alizations to other lattices, higher spins, etc. 

The nature of some of the various possible ground 
states of such a Hamiltonian are quite well understood. 
First, there are states that develop magnetic long range 
order and break the spin rotation symmetry. The sim- 
plest example (and the one that we will focus on) are 
coUinear antiferromagnets where the order parameter is 
a single vector (the Neel vector), defined to describe 
a state of staggered magnetization, 

Sr = £rNr, (1.2) 
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FIG. 1: Ground states of the square lattice S = 1/2 antifer- 
romagnet studied in this paper. The coupling g controls the 
strength of quantum spin fluctuations about a magnetically 
ordered state, and appears in Eq. 12. H (the classical limit 
is g — 0). There is broken spin rotation invariance in the 
Neel state for g < Qc, described by the order parameter Nr 
in Eq. 11.21 . The VBS ground state appears for g > gc, and 
is characterized by the order parameter tpyss in Eq. 11.41 - 
the distinct lines represent distinct values of {Sr ■ Sri) on each 
link. The VBS state on the left has "columnar" bond order, 
while that on the right has "plaquette" order. The theory 
in Eq. 11.711 applies only at the QCP g = g^ at its critical 
point obtained at s = Sc. 



where 



(—1)^"'"*' 



(1.3) 



is +1 on one checkerboard sublattice and —1 on the 
other. The Neel state has (Nr) ^ and independent 
of r (see Fig^, but more generally iV^ is presumed to 
vary "slowly" on the lattice scale over at least most of 
space. The low energy excitations of the Neel state are 
simply linear dispersing spin waves. 

It is now recognized that a variety of quantum param- 
agnetic ground states are also possible where quantum 
fluctuations have prevented the spins from developing 
magnetic long range order, and so {Sr) = 0. Such para- 
magnetic states can be broadly divided into two groups. 
First, there are states that can be described as 'Valence 
Bond Solid' (VBS) statesii. In a simple caricature of 
such a state, each spin forms a singlet with one particu- 
lar other spin resulting in an ordered pattern of 'valence 
bonds'. For spin-1/2 systems on a square lattice, such 
states necessarily break lattice translational symmetry. 
The so-called "columnar" and "plaquette" ordering pat- 
terns (see Fig. ^ are described by a complex VBS order 
parameter V'vbSi where 



Sr-Sr+i ^ RelV'VBsK-l)", 
Sr ■ Sr+y ^ Im[V'VBs](-l)^, 



(1.4) 



and r = (x, y) (here columnar states have V'vbs ^^^^ 
and positive, while plaquette states have V'vbs ^^^^ ^^'^ 
negative). In these states there is an energy gap for 
spin-carrying S = 1 quasi-particle excitations; these 



'triplons'— are quite distinct from spin waves, and are 
instead adiabatically connected to spin excitons in band 
insulators. A second class of more exotic paramagnetic 
states are also possiblei^iiiiiSiiSiil in principle: in these 
states the valence bond configurations resonate amongst 
each other and form a 'liquid'. The resulting state has 
been argued to possess excitations with fractional spin 
1/2 and interesting topological structure. 

Our focus will be on the nature of the evolution of the 
ground state between these various phases. Our primary 
example is that between the ordered magnet and a va- 
lence bond solid. We also discuss the phase transitions 
between valence bond solid and 'spin' liquid phases (see 
Section FVllip . Qualitatively similar phenomena will be 
ahown to obtain at both these transitions. 

Both the magnetic Neel state, and the valence bond 
solid are states of broken symmetry. The former breaks 
spin rotation symmetry, and the latter the symmetry of 
lattice translations. The order parameters N and "iAvBS 
associated with these two different broken symmetries 
are very different. A LGW picture of the evolution be- 
tween these two distinct ground states would be formu- 
lated in terms of an effective action that is a functional 
of N and V'vbs- Such a construction would suggest ei- 
ther a first order transition, or passage through an in- 
termediate phase which breaks both kinds of symmetry. 
(Actually, the general LGW analysis also allows an in- 
termediate 'disordered' state with neither order, but this 
possibility was excluded in early analyse^i; this exclu- 
sion was already an indication that LGW theory did not 
apply here.) A direct second order transition would be 
expected only by further fine-tuning to special multicrit- 
ical points. Our central thesis is that this expectation is 
wrong. A generic second order transition is possible be- 
tween these two phases with different broken symmetries. 
The resulting critical theory is however unusual and not 
naturally described in terms of the order parameter fields 
of cither phase. Instead, the natural description is in 
terms of spin-1/2 "spinon" or CP^ fields Za (a = 1,2 is 
a spinor index). The Neel order parameter is bilinear in 
the spinous: 



N . 



(1.5) 



Here a is the usual vector of Pauli matrices and multipli- 
cation of the spinor index is implied. The fields Zq, create 
single spin-1/2 quanta, "half" that of the spin-1 quanta 
created by the Neel field N . 

The spinon fields Za so defined have a U(l) "gauge" 
redundancy. Specifically the local phase rotation 



(1.6) 



leaves the Neel vector invariant and hence is a gauge 
degree of freedom. Here r is the imaginary time coordi- 
nate. Thus the spinous are coupled to a U(l) gauge field 
a^i(r, r) (we will use the Greek indices /i, i^, . . . to repre- 
sent the three spacetime indices x, y, r). Our central the- 
sis - substantiated by a variety of arguments to follow - 
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is that the critical field theory for the Neel-VBS transi- 
tion is just the simple continuum action — J cPrdrCz, 
and 

N 
a=l 

+ K{ef^^i^d^ai^f , (1.7) 

where = 2 is the number of z components (later we will 
consider the case of general N), |zp = J2a=i l-^aP, and 
the value of s is to be tuned to a critical value s = Sc so 
that Cz is at its scale- invariant critical point. The same 
action with a simple modification also describes the crit- 
ical field theory for systems with easy-plane anisotropy, 
with the addition of the simple term 

^ w\zi\^\z2\^ (1.8) 

with w < 0. We will discuss in more detail later why 
these would describe stable critical points - perhaps the 
most direct evidence comes from the numerical simula- 
tions reported in Ref. I21I of a lattice model of a CP^ field 
coupled to a noncompact gauge field (a lattice version of 
Eqn. II. 7|) . where a continuous transition was found in 
both the isotropic and easy plane cases. 

How can this action describe the onset of VBS order 
when it does not contain "^vbs 7 and the Za are closely re- 
lated to the Neel order parameter? In writing Eq. H1.7|l . 
we have tacitly assumed to be a single-valued contin- 
uous field. In a more careful lattice implementation of 
Eq. 1)1. 5|l . however, the resulting gauge field that appears 
is compact, i.e. defined only modulo 27r. This allows for 
the presence of topological defects occurring at a single 
instant of space-time ( "instantons" ) called monopoles, at 
which magnetic flux dxCLy — dyQx is created or destroyed 
in integer multiples of 2tt. In general, Eq. (|1.7|l should 
thus be supplemented by terms which create or destroy 
such 2-K fluxes, or equivalently insert monopoles into the 
partition function: 

00 

>C„.p = 5]A„(r)(KJ" + bt^n, (1.9) 

n=l 

where uj^. and v^.^ insert monopoles of strength 27r and 
— 27r at the space-time point (r, r), respectively. Remark- 
ably, it has been shown by Read and SachdeAiii*i^ that 
this operator may be identified with the VBS order pa- 
rameter, i.e. 

Vrr ^ 1pVBsir,T). (1-10) 

A simple argument to this effect will be given in Scc. lIIII 
Thus VBS physics is implicitly (albeit highly non- 
linearly) contained in the gauge theory of Eqs. 1)1. 7|) and 
(|1.9() . Our claim that Eq. (|1.7() without the monopole cre- 
ation terms of Eq. ^l.Up describes the critical properties 
of the Neel-VBS transition requires that the monopole 
"fugacities" X„ are irrelevant in the renormalization 



group sense at the QCP. Later sections of this pa- 
per will give a variety of compelling arguments, relying 
upon destructive quantum interference between different 
monopole events, for this irrelevance for spin-1/2 antifer- 
romagnets. The arguments are based on quantum Berry 
phase effects described first by HaldaneiS, which render 
Xn{r) oscillatory and negligible for n ^ (mod 4) for spin 
S = 1/2 (a different derivation appears in Appendix IXI 
and in the review in Ref. l20|) . 

Although monopoles can be neglected at the QCP, this 
is not true at low energies in the VBS phase. Indeed, it 
is well-known from studies of pure compact U(l) gauge 
theories, that the fugacities A„ are always relevant in 
the absence of gapless "matter fields" {i.e. the Za), so 
that monopoles inevitably proliferate in this case. This 
proliferation leads to a "condensation" of the monopole 
operator, (vrr) ^ (V'vbs) Oi hence VBS orde r^"^'^^ . At 
the same time it generates a gap for the gauge "photon" . 
In renormalization group terminology, the monopole con- 
densation in the VBS phase - despite the fact that the 
A„ are negligible at the QCP - indicates that (some) A„ 
are "dangerously irrelevant" . 

It is important to note that such monopoles have a 
natural topological interpretation in terms of the confor- 
mations of the Neel ordered state. In particular, low but 
non-zero energy configurations of the antiferromagnet are 
described by states with slowly-varying Neel vector (at 
least at spatial infinity) of constant amplitude, 

Nr = \N\nr. (1.11) 

Such classical configurations with finite energy admit 
topological defects known as skyrmions (see Fig|2l. The 
total skyrmion number associated with a configuration 
defines an integer topological quantum number Q: 

Q = ^ J (frn ■ dxfi X dyh. (1-12) 

Remarkably (see Sec. lIIII and Ref. 12^, the skyrmion den- 
sity is simply related to the magnetic flux of the gauge 
field a^, 

2tiQ = j (fx{dxay - dyax). (1-13) 

Thus the monopole instantons that change the gauge flux 
by ±27r describe events in which the skyrmion number 
changes by ±1. Thus the flux creation operator vl^ can 
also be interpreted as a skyrmion creation operator. The 
skyrmion number changing events may be represented 
graphically as "hedgehog" configurations of the Neel vec- 
tor in space-time (See Fig. O. The irrelevance of the 
monopole fugacities at the Neel-VBS critical point is thus 
equivalent to the irrelevance of hedgehog fugacities in a 
semiclassical description. Further, the recognition that 
such monopole events can be characterized as changing 
(as a function of time) the skyrmion number Q enables 
another interpretation of their irrelevance. In particular, 
in the critical fixed point theory in Eq. (|1.7|l - i.e. at 
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FIG. 2: A skyrmion configuration of the field hr- In (a) 
we show the vector {n^,n^) at difi'erent points in the XY 
plane; note that h oc [ — 1]^^^ Sr, and so the underlying spins 
have a rapid sublattice oscillation which is not shown. In (b) 
we show the vector {n^,n^) along a section of (a) on the x 
axis. Along any other section of (a), a picture similar to (b) 
pertains, as the former is invariant under rotations about the 
z axis. The skyrmion above has n(r = 0) = (0,0,1) and 
n(|r|^oo) = (0,0,-1). 



low energies near the QCP - the skyrmion number Q is 
strictly conserved. The emergence of this conserved topo- 
logical quantum number is the most fundamental mean- 
ing of the irrelevance of the instantons. 

We will also use this emergent topological conservation 
law as a definition of a "deconfined" QCP. Indeed, typi- 
cally the gauge theories that arise in various slave particle 
descriptions of quantum magnets are compact. Specializ- 
ing to a U(l) gauge theory, the compactness means that 
instanton or monopole events in which the magnetic flux 
changes by 2tt are allowed configurations of the gauge 
field in space-time. The proliferation of these instanton 
events leads to confinement of the slave particles in the 
gauge theory. In contrast, in a non-compact theory - 
which emerges at low energies when monopoles are ir- 
relevant - the total magnetic flux is strictly conserved. 
This is a topological conservation law and may be under- 
stood as a global U(l) symmetry in an appropriate dual 
description. Indeed, we will explicitly construct such a 
dual theory for the case of easy-plane anisotropy (and in 
some other related models). Quite generally, then, the 
emergence of a non-compact U(l) gauge theory at the 
critical point between the Neel and VBS phases signi- 






(b) 



FIG. 3: A monopole event, taken to occur at the origin of 
spacetime. An equal-time slice of spacetime at the tunnelling 
time is represented following the conventions of Fig|21 So (a) 
contains the vector {n^,ny); the spin configuration is radi- 
ally symmetric, and consequently a similar picture is obtained 
along any other plane passing through the origin. Similarly, 
(b) is the representation of (n^,n^) along the x axis, and a 
similar picture is obtained along any line in spacetime passing 
through the origin. The monopole above has hr = r/\r\. 



ties an extra emergent (dual) global U(l) symmetry for 
the critical theory that is not present in the microscopic 
Hamiltonian. This provides a rather precise characteri- 
zation of a 'deconfined' critical point. 

An important property of the deconfined fixed points 
discussed in this paper is the appearance of two distinct 
diverging length (or equivalently two time) scales close 
to the transition - one of which rises as a power of the 
other. This is directly due to the dangerous irrelevance 
of monopoles. For the Neel- VBS transition on approach- 
ing from the VBS side there is of course a diverging spin 
correlation length ^. However just beyond this length 
scale the system has not yet chosen to pin itself into any 
particular VBS ordered state. Rather it may be char- 
acterized as fluctuating between different VBS configu- 
rations. It settles down to a particular ordered state at 
a larger length scale ^vbs- This new length scale may 
also be characterized as the thickness of a domain wall in 
the VBS order. The universal crossovers associated with 
the critical fixed point describe the behavior on passing 
through the length scale ^. These are described by the 
critical theory in Eq. (|1.7|) . As explained above, this 
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critical theory is monopole-free. The second crossover 
associated with the length scale fvBS describes how the 
system evolves from the paramagnetic phase associated 
with the monopole-free theory Eq. 1)1. 7|l to the true VBS 
phase that obtains when monopoles eventually prolifer- 
ate. Further details of the physics at the scales ^ and fvBS 
appear in Section rVIII where we also show in Eq. (|7.3|l 
that ^vBS diverges as a power of ^ which is greater than 
unity. 

Over the last several years we have become familiar 
with the notion of fractionalization of quantum numbers 
in stable phases in condensed matter. In contrast, the 
fractionalization phenomena obtained in this paper are 
specific to the critical point separating two conventional 
phases. These 'fractional' particles - the spinous - are 
not present (i.e confined or condensed) at low energies 
on either side of the transition but appear naturally at 
the transition point. Likewise the emergent gauge field 
that mediates interactions between the fractional parti- 
cles is also specific to the critical point. On approaching 
the critical point, the confinement length scale diverges. 
Thus 'deconfincmcnt' appears right at the transition. 

We will also briefly discuss the phase transitions be- 
tween different quantum paramagnetic ground states. In 
particular we will argue that the existing theory for the 
transition between a VBS state and a fractionalized spin 
liquid implies that the corresponding critical point is also 
described by a deconfined U(l) gauge theory in pr ecisely 
the same manner as above. Furthermore, Refs. I23l24l 
argue that (at least under certain conditions) there are 
direct transitions between two different VBS phases that 
are also described by deconfined critical points with a 
U(l) gauge structure. 

There are several general lessons to be learnt from the 
results in this paper. First, we see that two dimen- 
sional spin- 1/2 quantum magnetism is full of examples 
of 'deconfined' quantum critical points which contradict 
the LGW paradigm for critical phenomena. This sug- 
gests that in more complex quantum systems {e.g. with 
fermions or disorder) novel critical phenomena may well 
be quite commonplace. Such deconfinement may be at 
the root of interesting non- Fermi liquid critical phenom- 
ena observed in the heavy fermion materials and possi- 
bly in the cuprates as well. Second, our results resolve 
some long-standing controversies in the field of two di- 
mensional quantum magnetism and have direct implica- 
tions for experiments and numerical work in the field. 

Third, our results shed some light on questions of con- 
finement in gauge theories in two spatial dimensions. It 
was shown by Polyakov several years age^ that in two 
spatial dimensions for pure gauge theories (i.e without 
any matter fields) instantons generically always prolif- 
erate and drive the theory into a confined phase. The 
behavior in the presence of dynamic matter fields (par- 
ticularly with fermionic matter) is much less understood 
and is a subject of some controvers}i2Si2L22i224ffi. The re- 
sults in this paper show that with bosonic matter there 
are at least isolated critical pointsSS at which deconfine- 



ment is obtained (and the instantons disappear at long 
scales). While typically reaching criticality in a bosonic 
system requires some fine-tuning, fermionic systems can 
have stable critical phases. This supports the specula- 
tion that stable deconfined phases exist in two dimen- 
sional compact U(l) gauge theories coupled to fermionic 
matter^^ . If true this would have interesting implications 
for the theory of spin liquid phases of quantum spin sys- 
tems. These points are discussed further in Section [IV Bl 
and Appendix IFl 

Apart from these general notions, there are also a num- 
ber of specific physical ramifications of the proposed crit- 
ical theory for the Neel-VBS transition. One immedi- 
ate consequence is that the anomalous dimension of the 
magnon operator is much larger than is usual at D — 2+1 
dimensional fixed points. Thus the magnon spectral func- 
tion will be extremely broad right at the critical point. 
Many other implications are explored in some detail later 
in this paper and summarized in the overview in Sec- 
tion Inl 



II. OVERVIEW 

In this section we provide an overview of the main ideas 
in this paper. 

A. History and Precedents 

We begin by recalling some important prior results in 
the theory of quantum magnetism on the two dimen- 
sional square lattice. In the Neel phase or close to it, 
the long distance low energy fluctuations (of the orien- 
tation) of the Neel order parameter are captured by the 
quantum 0(3) Non-Linear Sigma Model (NLctM) with 
the Euclidean action (we have promoted the lattice co- 
ordinate r = (x, y) to a continuum spatial co-ordinate, 
and r is imaginary time): 



So 
Sb 



So + Sb 



(2.1) 



Here oc e^S'r is a unit three component vector that rep- 
resents the Neel order parameter (the factor is defined 
in Eq. (|1.3|l ). The term Sb contains crucial quantum- 
mechanical Berry phase effects, and is sensitive to the 
precise quantized value, S of the microscopic spin on 
each lattice site: Ar is the area enclosed by the curve 
mapped out by the time evolution of firiT) on the unit 
sphere. These Berry phases play an unimportant role 
in the low energy properties of the Neel phase^^, but are 
crucial in correctly describing the quantum paramagnetic 
phas Q i Tir' T - . We will expand these earlier results to show 
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here that they also modify the quantum critical point 
between these phases, so that the critical exponents are 
distinct from the theory without Sb studied earhe u?^'?? . 

To understand the summation over r in Sb, recall 
that, as described in Section^ in two spatial dimensions, 
smooth configurations of the Neel vector admit skyrmion 
topological defects characterized by the integer topologi- 
cal charge Q. The Berry phase Sb vanishes^'^'' for all 
smooth configurations even if they contain skyrmions. 
For such smooth configurations, the total skyrmion num- 
ber Q is conserved. Crucially, however, the skyrmion 
number changing monopole events are not everywhere 
smooth. It was shown by HaldaneiS that the summa- 
tion over r in 5s is non- vanishing in the presence of 
such monopole events. Precise calculationi2*iS gives a 
Berry phase associated with each such skyrmion chang- 
ing process which, for 25* = 1 (mod 4), oscillates 
rapidly on four sublattices of the dual lattice (see Ap- 
pendix 0. This leads to destructive interference which 
effectively suppresses all monopole events unless they are 
quadruplediiiSiiS (i.e they change skyrmion number by 
four). 

The NLcrM field theory augmented by these Berry 
phase terms is, in principle, powerful enough to correctly 
describe the quantum paramagnet. Summing over the 
various monopole tunnelling events shows that in the 
paramagnetic phase the presence of the Berry phases 
leads to VBS orderiiiii^. This crucial result from prior 
work identifies the VBS phase as resulting from a prolif- 
eration of monopoles in the presence of non-trivial Berry 
phases. The non-trivial identification of the VBS order 
parameter expressed from bond energies in Eq. (|1.4I) as 
the skyrmion creation operator is remarkable. When this 
operator acquires an expectation value VBS order results. 
In this manner 5„ contains within it the ingredients de- 
scribing both the ordered phases of H. 

An important conclusion which follows from this body 
of work, via the above identification, is that a direct sec- 
ond order transition from the (coUinear) Neel phase to a 
translation symmetric 'spin-liquid' is likely to be absent 
in two spatial dimensions. This still leaves several pos- 
sibilities for the evolution of the ground state from the 
Neel to the VBS phases. For instance, there could be two 
transitions with an intermediate phase that breaks both 
Neel and lattice symmetries (an intermediate phase that 
breaks neither symmetry is excluded), or simply a first 
order transition. 

The possibility of a direct second order transition be- 
tween Neel and VBS phases is hinted at by several re- 
sults in the existing literature. First, note that in the 
Neel phase monopole tunneling events are absent at long 
length and time scales. In the quantum paramagnet these 
monopole tunneling events have proliferated. The Hal- 
dane phases then lead to VBS order. The existence of 
a monopole condensate is clearly incompatible with long 
range Neel order. Thus to the extent that the broken 
lattice symmetry of the VBS state is a consequence of 
the proliferation of monopoles it competes with the Neel 



state. A direct transition from Neel to VBS then becomes 
conceivabl9ii4i2i2&. 

A second hint comes from examining large-A^ stud- 
ies of SU(A^) quantum spin modelsiiii^. In the limit 
N — oo (and in a bosonic representation of the spins) 
there is a second order transition between the Neel phase 
and a quantum paramagnet. In this limit the param- 
agnet breaks no symmetries (including lattice symme- 
tries). Furthermore, it supports gapped spin-1/2 exci- 
tations (known as spinous). However, both features are 
known to be artifacts of the limit N = oo. Upon includ- 
ing finite N corrections, broken lattice symmetry appears 
(leading to a VBS phase). The spinous also feel a gauge 
force that leads to their confinement and disappear from 
the spectrum. What about the vicinity of the transition? 
To answer this, it is instructive to examine the various 
length scales in the problem in the paramagnetic state. 
First, there is the spin correlation length that diverges on 
approaching the transition. Note that deep in the para- 
magnetic phase this length stays constant when N ^ oo. 
Then, there is the length scale at which VBS order ap- 
pears. As there is no VBS order at = oo this length 
scale must diverge as A^ ^ oo - hence it must be much 
bigger than the spin correlation length in the \a,Tge-N 
limit. Finally, there is a third (somewhat loosely defined) 
length scale that may be thought of as the length scale as- 
sociated with spinon confinement. Clearly this scale also 
diverges as A^ — > oo and is much bigger than the spin 
correlation length. Indeed calculations of the VBS and 
confinement length scales in the large- A^ limit show that 
they are ©(f^^), where ^ is the spin correlation lengthS 
and c is a constant. 

This suggests the possibility of a direct transition be- 
tween Neel and VBS states in the large-A^ limit where 
the monopoles (and hence their Berry phases) are irrele- 
vant at the critical fixed point, but are important in pro- 
ducing VBS order and confinement in the paramagnetic 
state. In critical phenomena parlance, the monopoles are 
dangerously irrelevant at the critical fixed point. 

A picture similar to this was in fact proposed several 
years ago by Chubukov et alM: However, it was not 
appreciated that the quadrupling of the monopoles, in- 
duced by the Berry phases, renders both the monopoles 
and their Berry phases irrelevant at the critical point 
(the distinction between the relevance of single versus 
quadrupled monopoles is absent in the large A^ limit'^^). 
In particular, it was assumed that the confinement length 
scale will stay finite at the transition, which was then 
modelled (for physical SU{2) spins) by the LGW theory 
obtained simply by neglecting Sb in Eq. H2.1|l : this is 
the 0(3)-invariant Wilson-Fisher fixed pointS. In light of 
the discussion above, it is clear that as the confinement 
goes hand in hand with the VBS order both confinement 
and VBS length scales diverge at the transition. Thus we 
might expect 'deconfinement' to appear at the transition. 

A weakness in the arguments of Chubukov et al.^^ was 
pointed out by Sachdev and Park—. The latter authors 
argued that there was a finite density of monopoles in 
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space-time right at the critical point of the 0(3) LGW 
model, and the Berry phases then implied the presence 
of finite VBS order at any such critical point. Based 
on this they suggested that a possible evolution between 
the Neel and VBS phases was through a region of coex- 
istence of both broken symmetries. However, they left 
open the possibility of a direct second order transition 
between the Neel and VBS phases, but argued that any 
such transition could not be described by the 0(3) LGW 
model. 

Our discussion here makes it clear it is necessary that 
the corresponding fixed point have no monopoles at 
long scales. The natural candidate is then precisely the 
fixed point governing the transition in the model with 
monopoles forbidden. The arguments of Motrunich and 
VishwanathSi, and our present analysis, on such models 
show that the appropriate critical theory is that of Lz 
in Eq. H1.7|l . It must be kept in mind that this critical 
theory is entirely distinct, with all critical exponents dif- 
ferent, from the 0(3) LGW model obtained by dropping 
Sb from Eq. I|2.1|l . The first indication that such a dis- 
tinct continuous transition could exist in the monopole 
suppressed 0(3) NLctM was from the work of Kamal and 
Murthji22. Recently, the transition in this model with 
monopole suppression was studied in Ref. |^ where a 
new approach that sidestepped the potential problems of 
Ref. was used. A continuous, non-Heisenberg tran- 
sitions with properties consistent with those of Ref. 
was found. Moreover, an independent numerical simula- 
tion of a CP"'^ model with a noncompact gauge field was 
performed (essentially Eq. 11.71) which also yielded a con- 
tinuous transition and exponents consistent with the sim- 
ulations of the monopole suppressed 0(3) NLaM. This 
provided a nontrivial check of both the essential correct- 
ness of the numerical calculations and direct support for 
the identification of Eq. H1.7|) as the critical theory for 
the monopole suppressed 0(3) NLaM transition. The 
easy plane deformation of these models was also stud- 
ied in Ref. l2lL where again a continuous transition was 
obtained. This transition was argued to possess the re- 
markable property of being self- dual. 

The possibility of deconfinement of spinous at the crit- 
ical point between Neel and VBS phases is also hinted 
at by a different consideration that is again motivated 
by the large- TV calculations. The excitations of both the 
Neel and VBS phases are conventional {i.e do not con- 
tain any fractionalized spinons). In a Schwinger boson 
description in terms of spin- 1/2 spinons this is achieved 
through confinement. However the detailed mechanism 
of such spinon confinement is different in the two phases. 
In the Neel phase (described as a spinon condensate) con- 
finement is achieved through the usual Higgs mechanism. 
On the other hand, in the VBS phase confinement is 
achieved through proliferation of instantons. This dif- 
ference in the confinement physics then makes it con- 
ceivable that neither mechanism is actually operational 
at the critical point and deconfinement obtains. 



B. Numerics 

There have been a large number of numerical stud- 
ies of the destruction of Neel order in the S = 1/2 
square lattice antiferromagnet^. While there is evi- 
dence for the existence of VBS order in the paramagnetic 
phase(^i*^2iiM4il^ the nature of the transition between 
the Neel and paramagnetic states has been difficult to 
address. A major obstacle is the well-known 'sign' prob- 
lem, which prevents large-scale Monte Carlo simulations. 
Until recently, all large scale studies of the destruction of 
Neel order have been on models with an even number of 
5 = 1/2 spins per unit cell, with a paramagnetic phase 
which does not break any lattice symmetries4S41 

The first large-scale study of the destruction of Neel 
order in a, S = 1/2 square lattice antiferromagnet, in a 
Hamiltonian which preserves a single 5 = 1/2 spin per 
unit cell and the full symmetry of the square lattice, was 
that of Sandvik et almt. This was on a model with a strong 
easy-plane anisotropy. Such easy plane models have been 
studied analytically previousl}i2Sii£ and will be pursued 
further in the present paper. Ref. [Bi found convincing 
evidence for VBS order in the paramagnetic phase. Fur- 
thermore, the VBS and Neel order appear to vanish at 
points close to each other, suggesting a direct second or- 
der transition in the class discussed in the present paper. 

We also note the wavefunction Monte Carlo work of 
Capriotti et alm^ on the SU(2) 5=1/2 antiferromag- 
net on the square lattice with first and second neighbor 
exchange. They found a 'resonating valence bond' wave- 
function characteristic of a spin liquid state. Our results 
here suggest that they were perhaps observing the decon- 
fined state characteristic of the critical point, and that 
they had not yet reached the crossover to VBS order at 
the longest scales. 



C. Plan of Attack 

In this paper, the proposal of a deconfined continuous 
Neel- VBS transition (as well as a VBS-spin liquid transi- 
tion) is substantiated by a variety of arguments. First, in 
Section imi we consider a concrete lattice CP^~^ model 
which, for N — 2, embodies the physics of the Neel state, 
the monopoles and their Haldanc Berry's phases (focus- 
ing on 5 = 1/2), and the VBS state. This model, intro- 
duced by Sachdev and Jalaber1i2&, and referred to here 
as the SJ model, provides a convenient starting point 
for theoretical analysis of the SU (2) invariant critical re- 
gion. We address the nature of the physically interesting 
N = 2 case by showing that, in the two limits = 1 (Sec- 
tion llV A|| and N = oo (Section llVBIl ). this model indeed 
sustains a deconfined critical point in the precise sense 
defined above. For A^ = 1, this can be directly shown 
using lattice duality transformations, which demonstrate 
an exact equivalence of the SJ model to a D = 3 classi- 
cal XY model with a fourfold symmetry breaking term 
which corresponds physically to strength four monopoles. 
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Such four-fold anisotropy is known to be irrelevant at the 
D = 3 XY transition^", establishing the deconfinement 
of this case. For N — oo, the scaling dimension of the 
4-skyrmion creation operator was computed previously 
by Murthy and Sachdev"^^, and is such that monopoles 
are again irrelevant. Hence we expect by continuity that 
monopoles are irrelevant for all N, including the inter- 
esting case N = 2. 

Second, in Sectional we consider specifically TV = 2, in 
the presence of additional (strong) easy-plane anisotropy. 
In this case, the SJ model may be rewritten as a pair of 
0(2) rotors (the phases of Zq) interacting with a compact 
U(l) gauge field. The latter may be analyzed directly us- 
ing duality techniques f Appendix IB|) . We obtain in this 
way an explicit dual representation in terms of complex 
"vortex" annihilation operators ipa (o: = 1,2) and a dual 
non-compact gauge field v4^, whose flux represents the 
(exactly) conserved uniform spin density S^. One may 
understand the relation to the CP^ variables by recogniz- 
ing that ip\ creates a +2tt vortex in Z2, while ipl creates 
a — 27r vortex in zi, both of which create physical 2t: vor- 
ticity in n~ = — in^ = 2;*Z2- The dual theory, £duai 
for ipa and is presented in Eq. H5.6|l . 

The dual representation £duai has an appealing semi- 
classical interpretation, described in detail in Sec. IV Al 
Briefly, the two types of vortices correspond to "merons" 
(half-skyrmions), in which the Neel vector points cither 
up or down inside the vortex core. The skyrmion number 
changing monopole events thereby correspond precisely 
to an event in which a vortex core tunnels from the up 
to down staggered magnetization or vice- versa. 

The advantage of this representation is that the 
(quadrupled) monopole fugacity appears explicitly as a 
local operator in terms of the vortex fields. Remarkably, 
if this fugacity. A, is set to zero (as appropriate at the 
QCP provided it is, as we argue, irrelevant), the dual 
action in Eq. H5.6() has precisely the same form as the 
original one, Eq. H1.7|l . More precisely, an exact equiv- 
alence can be demonstrated between lattice regulariza- 
tions of the original and dual theory in the absence of 
monopoleaSi. Thus, as found in Ref. |^ the proposed 
critical theory in the easy-plane case has an unusual self- 
duality property. 

The irrelevance of monopoles can then be argued in 
several ways. First, using the self-duality, each power 
of the skyrmion creation operator has the same correla- 
tions at the QCP as the corresponding power of the XY 
staggered raising operator n+. At the deconfined criti- 
cal point, fluctuations of ?i+ are expected to be stronger 
than they are at a conventional (confined) XY critical 
point. The corresponding quadrupled operator is already 
irrelevant in the latter case (as mentioned above), so we 
expect the four-skyrmion fugacity to be only more irrele- 
vant around the deconfined critical theory. This expecta- 
tion is supported by an explicit calculation in a large- 
(different from the N in the SJ model) generalization of 
the dual critical theory in Section IVI Al Further argu- 
ments are given in Appendix lUl 



In Sect ion llX Al we demonstrate for the easy plane case 
a direct derivation of the dual critical theory from a mi- 
croscopic bosonic representation of the underlying XY 
model, without utilizing either the NLcrM or SJ models. 

In Section IVIIII we show that analogous deconfine- 
ment obtains for a VBS to spin-liquid transition. The 
latter has already been discussed by several authorsiL^, 
and shown to be equivalent to the transition in a fully 
frustrated quantum Ising model, which has a simple XY 
critical fixed point unaffected at low energies by an irrel- 
evant 8-fold symmetry breaking term. We show that this 
description is in fact dual to a deconfined gauge theory 
in the same sense as above, and that the (dangerously) 
irrelevant 8-fold symmetry breaking term can likewise be 
interpreted as an irrelevant monopole fugacity. 



D. Organization of paper 

We will begin in Section imi with a general discussion of 
the important symmetries of the Hamiltonian, and their 
action on a variety of order parameters and operators. 
This section will also introduce the SJ model. The so- 
lution of the SJ model in a variety of tractable limits 
appears in Sections IIVI and FV Bl Section also contains 
a general, semiclassical description of the physics in the 
easy plane limit. The nature of the second-order critical 
point between the Neel and VBS states is discussed in 
Section lvTl A variety of predictions for the critical prop- 
erties on the Neel- VBS transitions with SU(2) and XY 
symmetries follow from our analysis. These are elabo- 
rated in Sec. IVIII readers not interested in the detailed 
theoretical analysis may skip ahead to this section with- 
out significant loss of continuity. Section FVIIII describes 
the deconfined critical point between the VBS and spin- 
liquid phases, as noted above. Finally, Section IIXI con- 
tains a variety of extensions of the results in this pa- 
per. Section lTX Al shows that the easy-plane Neel to VBS 
transition can be reinterpreted as a superfluid-insulator 
transition in an interacting boson system; the insulator 
in this case contains a density wave in the amplitude of 
the bosons to reside on 'bond' states. This approach 
also provides an alternative derivation of the dual model 
Eq. H5.6|l . Section IIX Bl briefly discusses the extension 
of our results to antiferromagnets with S > 1/2, while 
section HX CI considers the case of the honeycomb lattice. 
Section IIXDI discusses possible extension of our results 
to systems with Ising anisotropy. 



III. REPRESENTATIONS AND SYMMETRIES 

In this section, we describe the representation of mag- 
nets with local tendencies to Neel order in the NLcrM and 
CP^ (SJ model) representations. We describe the action 
of the physical symmetries on the corresponding h and 
Za fields in each case 
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We start from the action 5„ in Eq. I|2.1|l . The all- 
important Berry phase term in Sb is defined on the un- 
derlying square lattice, and it is clear that lattice scale 
cancellations are important for the physics we are in- 
terested in. It is therefore useful to return to a lattice 
formulation to obtain 



Sn — Sq + Sb 




We have now rewritten the areas Ar in terms of A, which 
represents the vector potential of a magnetic monopolc 
with flux Att placed at the center of fi space at each lat- 
tice site. This lattice model is a faithful representation 
of the original quantum antiferromagnet so long as g is 
large. The continuum limit of Sq in this model is clearly 
just what appears in Eqn. 12.11 The representation of 
the Berry phase used here leads directly to Eq. IjAip in 
Appendix ^ 

The Berry phases are crucial for a correct description 
of quantum paramagnetic phases. As described in the 
previous sections, it was shown by HaldanaiS that the 
Berry phases are non-vanishing only in the presence of 
monopole events. The calculations in Refs. 1&.1^ give 
the total phase (for spin- 1/2 magnets that we consider - 
for a derivation see here. Appendix 

nexp(z|CnAQ„) . (3.2) 

n 

Here the monopole is associated with a plaquette of the 
original lattice (or equivalently with a site of the dual 
square lattice), which is labelled by the index n. The 
product is over all locations of monopoles, and AQn = 
±1 is the change in skyrmion number associated with the 
monopole. Note that the periodic boundary condition 
along the time direction requires that the net change in 
skyrmion number is zero so that A(5„ = 0. The fixed 
integer field is 0, 1, 2, 3 depending on whether the dual 
lattice coordinate is (even, even), (even, odd), (odd, even) 
or (odd, odd), so that the phase factor associated with 
each monopole is l,i,— 1,— i on these sublattices (see 
Fig m in Section irvT)l . 

The oscillating nature of the Berry phase factors on ad- 
jacent plaquette leads to destructive interference between 
different tunnelling paths for single monopoles. Indeed 
this interference effectively kills all monopole events un- 
less they are quadrupled (i.e change skyrmion number 
by four). Hence only such quadrupled monopole events 
need be including in the quantum statistical mechanical 
partition sum. 

We have already indicated the remarkable identifica- 
tion of the VBS order parameter defined in Eq. H1.4|) with 
the skyrmion annihilation operator, ipVBS ~ as shown 



in Ref. This provide the crucial confluence of the 
loss of antiferromagnetic order (and consequent prolifer- 
ation of monopole events) with the onset of VBS order, 
counter to conventional LGW wisdom. Because of the 
importance of this result, we give a simplified derivation 
of this relation here. 

It is important to recognize that the VBS order param- 
eter in Eq. H1.4|l is entirely defined by its transformation 
under the symmetries of the Hamiltonian. Any other 
field with the same symmetry properties as V'vbs will, 
on general scaling and renormalization group grounds, 
be proportional to V'vbs in the critical region. Thus to 
prove the identification of V'vbs with the skyrmion cre- 
ation operator v, it is sufficient to show that the latter 
transforms identically to V'vbs under all symmetry oper- 
ations. 

As a topological index, the skyrmion number is un- 
changed under smooth global SU (2) spin rotations, hence 
the skyrmion number changing operator is also an SU (2) 
scalar. Likewise, V'VBS; being defined through scalar 
bond operators, is SU{2) invariant. Let us consider the 
effect of lattice symmetry transformations on v. In the 
functional integral this operator is defined by insertion of 
a space-time monopole. It is easy to see that under 7r/2 
rotations in the counter-clockwise direction about a di- 
rect lattice site (which we denote i?^/2)j the Berry phase 
associated with the skyrmion creation event changes by 
giTrS^ Thus if we denote by the skyrmion creation op- 
erator and specialize to S = 1/2, we have 

K/2 ■■ ivK (3.3) 

The skyrmion creation operator is actually defined on a 
plaquette - for the time being, we will label the plaquette 
by the lattice site at the top right corner. 

Lattice translation operations T^^y corresponding to 
translations by one unit along x, y directions of the mi- 
croscopic spin model have somewhat more subtle effects. 
First, in the rotor representation of n,., these translations 
are represented as 

Tx : fir -fir+x (3.4) 

Ty . Tl^ ^ Tij-^y (3.5) 

The change in sign of fi is due to the staggering implicit 
in its definition. Now note that the skyrmion number 
Q is odd under fi — > —h. Consequently Tx,y convert a 
skyrmion creation operator to an antiskyrmion creation 
operator at the translated plaquette. Furthermore due to 
the difference in the Berry phase factors for monopoles 
centered on adjacent plaquettes, there is a phase factor 
that is introduced by the translation. Simple calcula- 
tion gives the following transformation properties for , 
specializing again to S' = 1/2: 

Tx-vl -ivr+x (3.6) 
Ty -.vl +ivr+y. (3.7) 

It is now clear that a paramagnetic state with a uni- 
form expectation value of breaks these lattice symme- 
tries. For instance, if (v^) = {v) ^ 0, then R^/2TTx,Ty 
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are all broken. This suggests a plaquette ordered state 
such as that shown in the lower right of Fig. ^ A straight- 
forward comparison shows that, up to an innocuous con- 
stant pre-factor, the lattice transformation properties of 
V are identical to those of V'vbs determined from a more 
mundane analysis of Eq. H1.4(l . In particular, 

V ~ e-"/VvBS, (3.8) 

properly reproduces all the transformation properties of 
the VBS order parameter. Thus we may indeed identify 
the skyrmion creation operator with the order parameter 
for the VBS order. 

We have already introduced in the introduction the 
CP'^ "spinon" fields to represent the Neel order parame- 
ter. These may be introduced on the lattice, 

hr — zlaz^. (3.9) 

To maintain the unit magnitude of n^, the constraint 

1 - Zip + \z2\'^ = 1 should be imposed upon the spinor z = 
z{r,T) = {Z1,Z2)- 

One can show^^ that the partition function of the con- 
tinuum NLcrM (with the action in Eq. (|2.1I neglecting 
Berry's phase terms) is exactly reproduced by the con- 
tinuum CP^ model with the action 

Scp = J drd^l {d,, - laf,) z\\ (3.10) 

Here enters mathematically as a Hubbard- 

Stratonovich field, and by considering its quadratic 
Euler-Lagrange equation, one can deduce the relation of 
the skyrmion number Q to the gauge flux of a^j given in 

Eq. arm . 

As discussed above, to incorporate VBS phases it is im- 
portant to correctly account for the Haldane Berry phases 
associated with these instantons. An appropriate model 
has been constructed by Sachdev and Jalabert'^^. The 
Euclidean action of the Sachdev- Jalabert (SJ) model is 

Ss.i ~ Sz + Sa + Sb 

Sz = -i^z*,e"-z,+^,a + c.c., (3.11) 

i 

Sa = y ^ (e^iyAAi^OA - 27rg^) (3.12) 
Sb = i^^CriA^g^. (3.13) 

n 

Here we have put the complex spinon fields Zia on the 
sites, i, of a cubic space-time lattice in dimensions D = 

2- 1-1 (now, n denotes the sites of the dual cubic lattice) , 
and they satisfy a unit length constraints l^mp = 1 
on each lattice site. The represent the compact U(l) 
gauge field, and are defined on the links of the space- 
time lattice. Note that the Zia are minimally coupled 
to the gauge field. The term Sa represents the gauge 
field kinetic energy. The quantity is an integer gauge 
flux that is defined on the links of the dual cubic lat- 
tice. Its divergence which enters the term Sb represents 



the number of monopoles on the sites of the dual lattice. 
Consequently, Eq. (|3.13|) is identical to the contribution 
in Eq. H3.2I) . and Sb provides the Haldane Berry phase 
factors that make the action appropriate for describing 
spin-1/2 antiferromagnets on the square lattice. The Neel 
ordered phase is a 'Higgs' phase where the Zi have con- 
densed, while the VBS phase is a 'confined' phase where 
the Berry phases have led to broken lattice symmetry. 

The action 5s,i is clearly closely related to the lattice 
action 5„ in Eq. (|3.1|) . after replacing n by z via Eq. (|1.5|l . 
However, the corresponding Sb terms in Eqs. H3.1|l and 
(13.13(1 do appear rather different - they are related by the 
Berry phase summation carried out by Haldanei^. Here, 
we establish the connection between these two forms 
of Berry phases in Appendix ^ further details on the 
derivation of Eq. ((3.13|l from the microscopic antiferro- 
magnet appear in Refs. I 20l36l38i 

If, as we will argue in the following, monopole events 
can indeed be neglected at low energies near the QCP, we 
can set = 0. Taking then a naive continuum limit of 
Eq. ((3.11I3.13|I gives precisely the proposed field theory 
of Eq. 1)1. 7() . We will, however, work directly with the 
lattice SJ model including monopoles in several of the 
sections to follow. 

As in any critical phenomenon, symmetry plays a key 
role in the discussion of the Neel-VBS transition. We 
therefore list here the various physical {i.e. non-gauge) 
symmetries of the problem and their action upon the Neel 
and spinon fields. The only continuous physical symme- 
try is spin-rotational invariance, either SU{2) or U(l) in 
the case of easy-plane anisotropy. Under such rotations, 
the Neel vector h and spinon field Zq, transform as global 
vectors and spinors, respectively. The remaining unitary 
symmetries are discrete operations of the space group of 
the square lattice, and can be composed from tt/2 rota- 
tions, translations, reflections, and inversions. As above, 
we denote 7r/2 clockwise rotations (around a direct lattice 
site) by i?7r/2 and unit translations in x and y hy T^, Ty 
respectively. Reflections x — )■ —x or y — > ~y around a 
lattice plane {i.e. leaving a row or column invariant) are 
denoted TZx,Tiy, and inversions about a site by /. Fi- 
nally, there is a non-unitary time-reversal operation T, 
which as usual takes microscopic spins Sj. — ^ — Sf. The 
transformation properties of the Neel, spinor, and gauge 
fields are given in Table 

Using these symmetry properties, one can determine 
the operators of the field theory corresponding to physi- 
cally interesting microscopic quantities in a spin model. 
Some of these are tabulated in the second column of 
Table ^1 of Section IV Bl and describe general situations 
which allow for easy plane anisotropy on the underlying 
magnet. (The SU{2) symmetric situation may be ob- 
tained as a special case). These include the easy-plane 
and hard-axis components of the Neel order parameter, 
Ar± = ± i]\[v^ and the VBS order parameter. 
Also included are the easy-plane and hard-axis compo- 
nents of the uniform magnetization M, all three being 
conserved for SU (2) symmetry but only the latter be- 
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Operation 


Coordinates 


Neel 


spinor 




jauge 


w 


Rti/2 


Xi > €ijXj 


invariant 


invariant 


Ui 








Xi ~* Xi + 1 


n — » — n 


• y t 

Za *''"q/3^0 






• 




Xi ^ Xi 


invariant 


invariant 


ai 


—* —ai 




/ 


r —r 


invariant 


invariant 




^ ^x/y 




T 




n —n 


• y t 


^x/y 


^ ^x/y 


• 



1/4 
3/4 1/2 







1/4 
• • 
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TABLE L Transformations of the Neel and spinor fields on the 
square lattice under the discrete symmetry generators. Here 
i = 1, 2 = a;, y is a spatial index, and eij — ia^j is the fully 
antisymmetric rank two tensor. Coordinate transformations 
of the arguments of the fields have been suppressed. 
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FIG. 4: Specification of the fixed field i9 = -C/4. The filled 
circles are the sites of the direct lattice, and i} resides on the 
sites of the dual lattice. 



ing conserved with easy-plane anisotropy. Two currents 
are also of interest: the spatial current jf of the con- 
served Ising (hard-axis) magnetization, and the vortic- 
ity three-current j"^ - the latter being meaningful in the 
presence of easy plane anisotropy. Finally, we may con- 
sider the CP^ gauge three-current — ^^i^xdj^ax, which 
in the continuum theory is identified with the topologi- 
cal current j^*; ^ \(-pLv\n ■ df^n x di^h. In a microscopic 
model, using the transformation properties of these oper- 
ators, one can construct a (rather complex) superposition 
of three-spin operators with these same transformation 
properties. Consider for instance, the time component 
Jq = eijOittj. On a square plaquette with central coordi- 
nate r, number the sites starting at the uppper-left corner 
of the plaquette and moving clockwise as 1, 2, 3, 4. Then 
one has 



Jo (r) - {-iy[Si • ^2 X 5^3 - 52 • ^3 X ^4 

-1-53 • 54 X 5i - 54 • 5i X 52] , (3.14) 

where the (—1)'^ takes opposite signs on the two sublat- 
tices of the dual lattice. 



SJ model at N 



Consider first N = 1 where Zi = e'"^* is simply a com- 



plex number of unit magnitude. Then 



5, = -2< 



cos (Aff) — a) 



(4.1) 



where sum is over the links £ of the cubic lattice. We 
indicate spacetime 3-vectors here in bold face, and the 
discrete lattice gradient by A. As discussed by SJ, this 

= 1 model displays a transition between a Higgs and a 
translation broken phase. The latter has a four-fold de- 
generate ground state due to lattice symmetry breaking. 
Simple symmetry arguments suggest a transition mod- 
elled by a Z4 clock model - as the four fold anisotropy is 
irrelevant at the D — 3 XY fixed point, this is in the 
3D XY universality class. SJ also provided numerical ev- 
idence supporting this expectation. As shown below, all 
of this is readily established by a duality transformation 
of the A^ = 1 model. 

To dualize the A^ = 1 SJ action we use a Villain rep- 
resentation of the Sz term in Eq. 14.11) : 



3"! 



ij ■ (Act) - a) 



(4.2) 



IV. SJ MODELS 



One useful generalization of the SJ model is to allow 
a — 1 . . . N above, so that z is an A^-component com- 
plex vector of unit magnitude. The SU{2) spin model 
corresponds to A^ = 2. It will be possible to analyse the 
limits A^ = 1 and N = 00. As argued in Section FVII CI 
the N = 1 model actually may be realized in a spin-1/2 
model in a staggered Zeeman field. The large-A^ limit 
describes ordering transitions of certain SU(A") quantum 
antiferromagnets and is less directly physical. Its main 
utility is its tractability. Similar behavior in both ex- 
treme limits - in particular the irrelevance of monopoles 
in both cases - suggests the same is true for the models 
with intermediate A^. 



The integer valued field j represents the current of the z 
field. We also decouple the Sa term in Eqn. 13.121 by a 
Hubbard-Stratanovich field b to write 



2K' 



ib • (A X a — 27rq) 



(4.3) 



Here, and below, the leading sum in the action extends 
over all sites/links/plaquettes over the cubic lattice, as 
needed. Performing the sum over the integer field q, we 

get 



b - A?9 = B 



(4.4) 



with i?r = (see Fig 0} and B an integer. If we 

now integrate over (f), we get the current conservation 
condition 



0. 



(4.5) 



13 



This may be solved by writing 

j = A X A 



(4.6) 



with A an integer. Integrating over the gauge field a, we 
obtain 



A X B =j. 

This may be solved by writing 

B = A + Ax 
with X integer. The action then reads 



i(A X A)^ 

2t ^ ' 



^(A(x + ^)+A)2 



(4.7) 
(4.8) 

(4.9) 



The hard integer constraints on A, x may be softened by 
adding terms 



cos(27rA) — ^ A„ cos(27rnx) 



(4.10) 



We may now shift X ^ X = X + A-^A' = A + Ax. 
The A' field is massive and may be integrated out. The 
remaining action for the x reads 



S 



E 



-i cos(27rAx) 



\n cos[27rn(x — '&)] 



(4.11) 

This describes an XY model with various n-fold 
anisotropy terms of strengths An, The shift by d leads 
to rapid spatial oscillations of these anisotropy terms un- 
less n = (mod 4) . Near the critical point in the contin- 
uum limit, the leading non- vanishing anisotropy term is 
at n = 4. The critical properties are therefore that of an 
XY model with four-fold anisotropy A4. The latter has a 
scaling dimension A4 > 3, which renders it irrelevant at 
the D = 3 XY critical point. 5° 

An overly cautious reader may object that uncon- 
trolled approximations have been made in softening the 
integer constraints on the A, x fields. However, all ma- 
nipulations up to Eq. H4.9|) are exact, and from this 
point an exact world-line representation may be obtained 
by implementing the integer constraints using the Pois- 
son resummation formula. The latter representation 
clearly describes charged relativistic particles for which 
charge non-conservation events oscillate spatially unless 
the charge is changed in multiples of four. On univer- 
sality grounds, one expects this model to be in the same 
universality class as an XY model with 4-fold anisotropy. 
Though we will not pursue it, a similar exact duality can 
be performed on the iV = 2 SJ model in a world-line 
representation, and may be used to somewhat more rig- 
orously argue for self-duality of the critical theory in this 
case. 

The results above can be interpreted physically as fol- 
lows. Let us first consider the vortices in the z conden- 
sate. These will carry gauge flux that is quantized in 



units of 27r. Such a 27r flux can end at a space-time 
monopole. Hence monopoles act as sources of the vor- 
tices of the z field. The Berry phases imply that these 
monopole events are quadrupled so that only processes 
where four vortices disappear (or are created) together 
are important in the continuum limit. Now if we forbid 
monopoles by hand, then the usual duality arguments 
map the model to a global XY model in terms of the 
vortex fields. The dual global U(l) symmetry of this XY 
model is precisely associated with conservation of vor- 
ticity. Including monopoles (which act as sources for 4 
vortices) introduces a four-fold anisotropy on this global 
XY model. Such an anisotropy is irrelevant at the 31? XY 
critical fixed point. Thus monopole events are again irrel- 
evant and (in the original representation) a theory where 
the z boson is coupled to a non-compact U(l) gauge field 
describes the transition. 



B. SJ model at large 

Now let us consider N large. In the limit N ^ 00 
the gauge field is non-fiuctuating and can be taken as a 
classical 'background' in which the z particles move. The 
minimum energy saddle point corresponds to a^^ = (up 
to gauge rotation). The z bosons are gapped and free in 
the paramagnetic state, while they are condensed in the 
ordered state. Now consider the nature of both states, 
and the transition, upon including fluctuations in a 1/iV 
expansion. It is useful to discuss the effects of instan- 
tons separately from other fluctuations. Ignoring instan- 
tons, the 1/A'' expansion proceeds along standard lines. 
In the ordered state, the gauge fields acquire a mass by 
the usual Anderson-Higgs mechanism. The gauge fiux is 
quantized in units of 27r - the associated point defects are 
the large- avatars of the skyrmion described previously 
directly at = 2. However, on the paramagnetic side 
the gauge fields are gapless and describe a 'photon' which 
disperses linearly at low energies. The transition is de- 
scribed by a field theory of z bosons coupled minimally 
to a non-compact U(l) gauge field. This transition is 
second order with critical exponents that evolve contin- 
uously from their values aX N — 00. In particular, con- 
sider the gauge invariant physical spin operator (which 
is the appropriate generalization to large- A'^ of the famil- 
iar Neel order parameter at A^ = 2). This is bilinear in 
the z fields. At A^ = 00, the spin operator therefore has 
a large anomoulous dimension ry = 1. This will acquire 
(calculable) corrections^^^ of 0{1/N) upon considering 
finite but large A^. Hence ry will be large for large but 
finite N. 

Now consider including instantons. It is important to 
realize that the entire gauge action is of order A^ in this 
theory. Consequently)^ the 'bare' instanton core action, 
obtained by integrating out the z fields in the presence of 
a background instanton configuration of the gauge fields 
in space-time, is of order A^. Thus the bare instanton 
fugacity is small (exponentially small in A^). In the or- 
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dered state, the inclusion of instanton events means that 
point defects with quantized 27r flux are no longer stable. 
The physics in the paramagnetic state is more interesting. 
Here the instantons proliferate and lead to confinement 
of the gapped z bosons. Furthermore, the gaplcss photon 
(present in the non-compact model) is rendered unstable. 
The Haldane Berry phases associated with the instan- 
tons lead to lattice symmetry breaking. As explained 
in Section Uni this follows from the observation that the 
instanton operators transform non-trivially under lattice 
symmetries. Hence if they acquire an expectation value, 
lattice symmetry is broken. 

Now let us consider the effect of instantons at the tran- 
sition. ^From the discussion in preceding sections, it 
is clear that the crucial question is whether the four- 
monopole event is relevant/irrelevant at the fixed point 
of this non-compact model. The scaling dimension of 
the p-monopole operator in this model was computed 
by Murthy and Sachdev'^^. For p = 4, their results 
give a scaling dimension oc N . Hence the instantons are 
strongly irrelevanliSi for large N . 

This then implies that the critical point of the non- 
compact theory is stable to inclusion of instanton events, 
even though the states on both sides of the critical point 
are qualitatively changed. In particular, consider ap- 
proaching the transition from the paramagnetic side. 
The proliferation of instantons in the paramagnetic state 
had two effects - to confine the spinous and to pro- 
duce VBS order. The irrelevance of the instantons at 
the critical fixed point implies that both the VBS order 
and the spinon confinement disappear at the transition. 
We note that as the bare instanton fugacity is exponen- 
tially small in iV, this perturbative analysis of their rel- 
evance/irrelevance is sufficient to determine the nature 
of the transition. In particular the alternate possibility 
that there is a coexistence region with width shrinking 
to zero as — > oo appears unlikely at large- A?^. 

It is also useful to interpret the results above in the 
context of other recent discussions of instantons in the 
literature22iS2i^2i^. The strategy of these, and other 
works, is to integrate out the z bosons, and to work with 
an effective action for the gauge field. This action will be 
of order A^. Consequently, it seems reasonable to assume 
that the gauge field dynamics is described to leading or- 
der in l/A?^ by a Gaussian action (this is equivalent to the 
RPA approximation), and to address issues of instanton 
physics within this Gaussian gauge action. Such an ap- 
proach will correctly describe the qualitative physics of 
the paramagnetic state. For the critical point itself, the 
form of the Gaussian gauge action is determined by scal- 
ing to be22iSa (see also Eq. 

^Arao|X||aT(K)|^ (4.12) 

where K is the Euclidean 3-momentum, refers to the 
transverse part of the gauge field, and CTq is a universal 
constant associated with the universal critical conduc- 
tivity of the z bosons at the transition at A^ = oo. Note 



that this action is more singular than the usual 'Maxwell' 
action - this originates in the integration over the mass- 
less critical modes of the Za fields. The action for a single 
instanton can be calculated within this Gaussian approx- 
imation, and is of order A^ In i where L is the system size. 
This suggests that instanton-anti-instanton pairs inter- 
act logarithmically with each other. It also suggests that 
the effect of instantons could be captured by studying 
a classical three dimensional Coulomb gas of instantons 
with pairwise logarithmic interactions. If this gas is in 
a plasma phase, free instantons have proliferated. On 
the other hand, one might also conceive a different phase 
where instanton-anti-instanton pairs are strongly bound 
to each other. For the classical 3D Coulomb gas, exami- 
nation of this issu o^^i?? has led to the conclusion that the 
logarithmic interaction is screened at long length scales 
into a short-ranged interaction, by bound instanton pairs 
at shorter scales (however, it was noted^*^ that the 
screening could fail at fine-tuned critical points). This 
screening then forces proliferation of free instantons, so 
that the Coulomb gas is in a plasma phase. How are we 
to reconcile this apparently general conclusion with our 
claim that the instantons are suppressed at the Neel-VBS 
critical point? 

This conundrum is resolved as follows. The Gaussian 
'RPA' action does not properly account for the effects of 
highly non-linear perurbations such as instantons. This 
is already clear from the results of Ref. I^S Within the 
Gaussian theory, the action of a strength p instanton 
scales with p as p^. This would imply that the scaling 
dimension of the p instanton operator scale as p^ - this 
disagrees with results of Ref. '39, which obtained a highly 
non-trivial dependence on p. In other words, even in the 
large- A^ limit, the Gaussian action is not sufficient to cor- 
rectly calculate the scaling dimension of the instantons: 
the non-linear terms in the gauge action all contribute 
in determining the instanton action^ even N — co. 
More significantly, we can likewise conclude that the RPA 
treatment of instanton interactions by a simple pairwise 
interaction is inadequate. The true instanton gas (even 
in the large-A^ limit) has a rather specific structure of 
higher order interactions, some of whose features are uni- 
versally determined by the fact that they arose from in- 
tegrating out particular gapless critical modes. If we at- 
tempt to compute the screening of instanton interactions 
by integrating out bound instanton-anti-instanton pairs, 
effects which renormalize the screening length are intri- 
cately intertwined with those that shift the position of 
the critical point between the magnetic and paramag- 
netic phases. Indeed, fine-tuning to be at the critical 
point between the magnetic and paramagnetic phases is 
all that is needed to also suppress the instanton plasma 
phase, and the 'naive' conclusion that the instantons are 
irrelevant at this critical point'^'^ is correct. Appendix IfI 
considers a specific toy model for which these arguments 
can be demonstrated explicitly. 

As we noted at the end of Section HJ the above rea- 
soning may also apply to fermionic models which have a 
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line of critical points^SiSi: in this case, the suppression of 
instantons may occur along the entire line, and not just 
at an isolated point. 



V. SPIN MODELS WITH EASY PLANE 
ANISOTROPY 

An alternate and particularly fruitful deformation of 
the model is provided by the situation where there is 
some easy-plane anisotropy on the underlying SU(2) spin 
model. Such an anisotropy tends to orient the spins pref- 
erentially perpendicular to the z-axis in spin space. In- 
deed precisely such an easy plane spin-1/2 model with 
both two particle and four-particle ring exchanges has 
recently been studied numerically^. A direct transition 
between Neel and valence bond solid phases was found. 

Consider first the fate of the global symmetries in the 
presence of easy plane anisotropy. A U(l) subgroup of 
symmetry of spin rotations about the z-axis of spin still 
survives. In addition there are a number of discrete sym- 
metries. Either under a unit translation or time reversal 
(see Table nj, the Neel vector changes sign 

fir ~* —fir- (5.1) 

This may be combined with a U(l) spin rotation in the 
XY-plane which restores the sign of n^,n^ to simply 
change the sign of alone. Thus — > — is a discrete 
symmetry in the easy-plane case. 

Easy plane anisotropy is readily incorporated into the 
non-linear sigma model description in Eqs. (|2.1(l or H3.1(l 
as a term 



(5.2) 



with w < (this is clearly related to Eq. HLHIl ). The 
global U(l) symmetry simply corresponds to a uniform 
rotation of all the fi vectors about the z-axis. 
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FIG. 5: The 'meron' vortices in the easy plane case. There are 
two such vortices, ^1,2, and t/ji is represented in (a) and (b), 
while i/)2 is represented by (a) and (c), following the conven- 
tions of Fig 121 The i/'i meron above has n(r = 0) = (0,0, 1) 
and n(|r| — > co) — {x,y,0)/\r\; the il}2 meron has n{r = 0) = 
(0, 0, —1) and the same limit as \r\ — > 00. Each meron above 
is 'half the skyrmion in Fig|5| this is evident from a com- 
parison of (b) and (c) above with Fig|^. Similarly, one can 
observe that a composite of i/'i and ^p2 makes one skyrmion. 



A. Semi-classical analysis 

Let us first think classically about this easy plane 
model. By classical we mean to focus on time inde- 
pendent configurations of the n-field and to ignore the 
Berry phase effects. The classical ground state simply 
consists of letting n be independent of position and lie 
entirely in the spin XY plane. Topological defects in this 
ground state will play an important role. With the easy 
plane anisotropy, these are simply vortices in the field 
n"*" = Tlx + iriy. More precisely, on going around a large 
loop containing a vortex the phase of n+ winds around 
by 2^771 with m an integer. 

What is the nature of the core of these vortices? In 
the core the XY order will be suppressed and the n vec- 
tor will point along the ±z direction. In terms of the 
microscopic spin model, this corresponds to a non-zero 



staggered magnetization of the z component of the spin 
in the core region. Thus at the classical level there are 
two kinds of vortices depending on the direction of the n 
vector at the core (see Fig. [S)). Note that either kind of 
vortex breaks the Ising-like — > —71^ symmetry at the 
core. 

Clearly this breaking of the Ising symmetry is an ar- 
tifact of the classical limit - once quantum effects are 
included, the two broken symmetry cores will be able 
to tunnel into each other and there will be no true bro- 
ken Ising symmetry in the core. This tunneling is often 
called an 'instanton' process that connects two classically 
degenerate states. 

Interestingly such an instanton event is physically the 
easy-plane avatar of the space-time monopole described 
above for the fully isotropic model. This may be seen pic- 
torially. Pictorially each classical vortex really represents 
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half of a skyrmion configuration. Such half-skyrmions are 
known as merons. As shown in Fig.jS] the tunnehng pro- 
cess between the two merons is equivalent to creating a 
full skyrmion. This is precisely the monopole event. 

Now the results of Haldane imply once again that (in 
the continuum) monopole events are quadrupled. Thus 
the only tunneling processes that survive in the contin- 
uum limit are those in which four merons with core spins 
along one direction come together and collectively flip the 
orientation of their core spins to produce four merons of 
the opposite kind. 

B. Easy plane in the CP^ representation 

It is extremely useful to also consider easy plane 
anisotropy in the framework of the CP^ representation. 
In this representation, the easy plane anisotropy was 
already presented in Eq. (|1.8|l . Let us first translate 
the above classical thinking into the CP^ representation. 
Suppose for this purpose w is negative but small, so that 
states with \zi\ = \z2\ are favored but not rigidly en- 
forced. Clearly, the preferred uniform classical ordered 
states satisfy 

l(^i>|- 1(^2)1 y^O, (5.3) 

so that n"*" = z*Z2 is ordered and there is no average value 
of = l^ip — |z2p- Now consider vortex configurations. 
Clearly a full 2t: vortex in ?t.+ can be achieved by either 
having a 27r vortex in zi and not in Z2 or a 27r anti- vortex 
in Z2 and no vortex in zi . Far from the vortex core both 
fields will have equal amplitude, but in the first choice 
the amplitude of the zi condensate will be suppressed 
at the core but (^2) will be unaffected. Consequently 
= I Zip — |2;2p will be non-zero and negative in the 
core. The other choice also leads to non-zero which 
will now be positive. Thus we may identify the two kinds 
of meron vortices with 27r vortices in the spinon fields zi^2 
respectively. 

To explore this analytically, we consider the behav- 
ior now deep in the easy plane limit, in which = 
l^iip — \zi2\'^ ~ 0. Together with the CP^ constraint 
l-^iiP + 1-2221^ — 1, this implies fixed magnitude for each 
component of Zi, so we may write 

-^=6'"^-, (5.4) 

where (j)ia G [0, 2tt) is the phase of the spinon field. The 
"kinetic" term of the SJ model action in Eq. H3.11|l is 
then modified to 

5z = -i^cos(A(/)Q - a) , (5.5) 

e,a 

with the other terms (5a, 5b) given as before in 
Eqs. and |TT^ . 

It is very useful both for further insight and for con- 
crete calculations to explore a 'dual' representation which 



focuses on the meron vortex fields. Although the form 
of the dual action is dictated completely by the general 
considerations of the previous subsection, we provide an 
explicit derivation in Appendix ^ by proceeding as in 
Section IIV Al with the duality transformation. We ob- 
tain the dual action iSduai — J (Pfd'TCdued with 

0=1,2 

+Wd\i^i\^\H^ + - XRe[{rM\ (5-6) 

We have used the short-hand notation = + 
|'!/'2p- Here ^1^2 denote quantum fields that destroy 
meron vortices whose core points in the up direction for 
^pl and down for ip2. 

We now show how the dual action of Eq. H5.6|l can 
be understood entirely on general grounds. As usual in 
dual theories, the net vorticity is conserved, correspond- 
ing to the overall U(l) symmetry of Eq. H5.6|l . This 
symmetry is gauged by the non-compact vector poten- 
tial A^, as usual in dual descriptions of two dimensional 
bosonic systems. Physically, the gauge field is required 
to embody spin (S*^) conservation of the original model, 
jfi = f^nvxduAx/iT being the 3-current of . Hence the 
dual magnetic and electric fields correspond to the spin 
density and spin current, respectively. Minimal coupling 
of the vortex fields to also gives them proper loga- 
rithmic interactions and magnus force dynamics. 

Clearly under the discrete Ising-like — symme- 

try, the two vortices get interchanged, i.e tpi — > ip2 and 
vice-versa. The dual action must therefore be invariant 
under interchange of 1 and 2 labels. 

Finally, if monopole events were to be completely ig- 
nored {i.e disallowed by hand) the total skyrmion number 
must be conserved. As is apparent from Fig [SJ a com- 
posite of a ■01 vortex and a ^'2 antivortex is precisely a 
skyrmion configuration of the fi field. Thus we may view 
skyrmion number conservation as the conservation of the 
difference of the total number of either species of vortices. 
This implies the global U(l) symmetry 

V'l V'iexp(ie) 

V'2 V'2exp(-?gi) (5.7) 

where p is a constant. 

As discussed at length above, monopole events destroy 
the conservation of skyrmion number, and hence this 
dual global U(l) symmetry. However as the monopoles 
are effectively quadrupled due to the Berry phase terms, 
skyrmion number is still conserved (mod 4). Thus the 
dual global U(l) symmetry must be broken down to Z4. 

The dual Lagrangian £duai in Eq. (|5.6|) is the sim- 
plest one that is consistent with all these requirements. 
In particular, we note that at A = the dual global U(l) 
transformation in Eq. H5.7|l leaves the Lagrangian invari- 
ant. The A term breaks this down to as required. 
Thus we may identify A as the fugacity of the (quadru- 
pled) monopole tunneling events, A ~ A4 in Eq. (|1.9|) . 
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Actually this action was first derived by completely 
different means in Refs. lisi andlssl The discussion above 
is however more directly physical, and gives an interpre- 
tation of the A term and of the other symmetries of this 
dual action. 

An important coupling constant in the above dual ac- 
tion is Wd, which appears as a sort of "anisotropy" . The 
exact lattice duality in the appendix in fact leads to a 
"hard spin" (rotor) model in which j-^il = \ip2\ = 1- The 
above continuum theory is arrived by "softening" this 
constraint. However, it is clear that the appropriate sign 
of Wd to connect microscopically to the original SJ model 
is < (and large). While the above symmetry argu- 
ments do not specify this sign, the model with Wd > 
presumably corresponds to rather different physics, and 
has no clear connection to the original SJ model. See the 
discussion in Sec. lIXDl for a possible physical application 
of this case. 

Note that apart from the A term, Eq. IjS.Gfl has exactly 
the same structure as the continuum CP^ theory Cz in 
Eq. (|1.7|l in the presence of the easy plane anisotropy 
Eqn. H1.8() . As we will discuss below, the A term (which 
represents the quadrupled monopole tunneling events) 
are irrelevant at the QCP: consequently, the QCP has 
a self-dual structure. 

We should further note that there is no connection be- 
tween Wd in the vortex action and the analogous parame- 
ter w in the continuum limit of the SJ model. The latter 
is clearly simply related to the physical spin anisotropy, 
corresponding for w < to easy-plane and for w > 
to easy-axis anisotropy. The point w = describes the 
SU(2) invariant magnetic QCP. This does not correspond 
to = in the dual theory. Indeed, there may be no 
dual theory whatsoever for any but the easy plane case 
(though see Sec. IIX D|) . 

A list of the representation of physical operators of in- 
terest in the original and dual representations is given in 
Tableim In the dual vortex theory, the XY ordered phase 
is simply characterized as a dual 'paramagnet' where 
both the ■01^2 fields are gapped. On the other hand, true 
spin paramagnetic phases correspond to condensates of 
the fields V'1,27 which break the dual gauge symmetry. In 
particular if both ■0i and ■02 condense with equal ampli- 
tude (-01) — (V'2) 7^ 0, then a paramagnetic phase where 
the global Ising symmetry is preserved results. Note the 
strong similarity between the description of the phases in 
this dual theory with that in terms of the spinon fields of 
the CP^ representation if we interchange the role of the 
XY ordered and paramagnetic phases. This is a symp- 
tom of an exact duality between the two descriptions that 
obtains close to the transition. At this point, the two 
descriptions do not appear wholly identical, due to the 
Z4 symmetry-breaking term A not present in the CP^ 
theory, and the compactness of the CP^ gauge field not 
present in the dual one. As argued above, the two dif- 
ferences represent one and the same physics, since the 
vortex tunneling events generated by A represent exactly 
the non-oscillating four-monopole events allowed in the 
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TABLE II: Operators in the easy plane CP^ (column 2) 
and dual (column 3) representations corresponding to phys- 
ical operators (column 1), in the notations of Sec. IIIII Here 
we have introduced the classical gauge field configuration for 
a unit point fiux, with eijdiAj{x) = 2-k5'^(x). The symbol 
Ej represents the j- componenet of the electric field opera- 
tor that corresponds to the gauge field in the CP^ repre- 
sentation. We have also used the symbol 9^, defined by 

fdt^g = kifdi^g - {diif)g]- The symbol jf is the current 
of conserved magnetization, while and are the three- 
currents of vorticity and gauge flux, respectively. 

SJ model. In the next section, we will argue that these 
events are irrelevant in the scaling limit near the QCP, 
making the duality between the two descriptions com- 
plete. 

As indicated in Table ^1 the combination V'vbs ~ 
serves as the order parameter for the translational 
symmetry broken VBS gro und st ate. This may be seen 
from the analysis of Refs. '38^48?. Alternately this may 
be seen by the identification described in Sec. IIIII of the 
skyrmion creation operator with the order parameter for 
translation symmetry breaking. Such a condensate of 
"01,2 breaks the global Z4 symmetry of the action in 
Eq. H5.6|l . The preferred phase of V'vbs depends on the 
sign of A, the two inequivalent sets of preferred direc- 
tions corresponding to columnar and plaquette patterns 
of translational symmetry breaking. 

VI. PHASE TRANSITIONS 

A. Easy plane limit 

Consider the dual vortex action in Eq. IjS.Gfl . In mean 
field theory the transition happens when the parameter 
Sd becomes smaller than zero and can clearly be second 
order. Fluctuation effects will modify the mean field be- 
havior in important ways. Consider first the properties 
of the transition when A = 0, i.e in the absence of in- 
stanton events. The resulting model has recently been 
studied in Ref. |2l|. Remarkably, as argued there, the 
model has the property of being self-dual - the ordered 
and paramagnetic phases get interchanged under the du- 
ality transformation. To understand this first note that 
in the A = limit, the dual action Eq. H5.6|l has pre- 
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cisely the same structure as an easy-plane CP"'^ model 
with a non-compact U(l) gauge field (as in Eq. (ll.Tfl ). 
As this same limit actually corresponds to disallowing all 
monopole events, in the spinon description we must work 
with a non-compact gauge field. Then the exact same 
field theory obtains both in terms of the spinon fields z 
(in Eq. (^3) and in terms of the meron vortices ip in the 
easy plane limit (in Eq. H5.6|l 'l and ignoring instantons. 

It was established in Ref.[^via numerical Monte Carlo 
simulations that a continuous ordering transition exists 
in this model with the non compact gauge field. The fixed 
point controlling this transition in this limit is therefore 
described by a self-dual field theory. Note that in either 
representation the natural fields of the theory are not 
those associated with the 'physical' boson operator (ei- 
ther n"*" or the skyrmion creation operator). Rather the 
theory is expressed most simply in terms of 'fractional- 
ized' fields - namely the spinous or the meron vortices. 
In particular, the physical n"*" field is a composite of two 
spinon fields and likewise the skyrmion field is a compos- 
ite of the two meron fields. 

Let us now imagine including instanton events. This is 
most easily accessed in the vortex representation where 
it simply amounts to letting A 7^ 0. This is the main ad- 
vantage of the dual representation - the non-trivial non- 
local effect of instantons is represented as a simple local 
perturbation in the dual theory. We may now address 
the question of relevance/irrelevance of instantons at the 
A = fixed point. This is determined by the scaling 
dimension A of the (V'iV'2)^ = V'vbs operator, in prin- 
ciple determined from the two-point correlation function 
of this operator in the (non-trivial) theory with A = 0: 



non-compact gauge theory with Lagrangian 



V'VBS{^)V'VBS(^') 



1 



A=0 



(6.1) 



where x, x' are space-time coordinates. Hence A is deter- 
mined by the correlations of the fourth power of the phys- 
ical VBS order parameter, and one requires A > = 3 
for irrelevance. Being self-dual, the same anomalous di- 
mension should be ascribed to the correlations of the 
physical boson (XY order parameter). The A = criti- 
cal fixed point describes an XY ordering transition where 
the physical boson field is a composite of the fundamen- 
tal fields of the theory. We therefore expect that corre- 
lators of the physical boson (and its various powers) will 
decay with an anomalous dimension that is larger than 
the corresponding one for the ordinary XY transition in 
D = 2-\-l dimensions. Now for the usual XY fixed point 
four-fold symmetry breaking perturbations are known to 
be irrelevant, i.e. have a scaling dimension A4 > 3. This 
then implies that a small A will be irrelevant by power 
counting at the A = fixed point of the present model as 
well (see also Appendix Q. 

This latter expectation can be checked in an appropri- 
ate large N generalization. In particular, consider the 



2N 
i=\ 

N 

+ tiNie^^^xdM^ - A ^[(7A*^.+Ar)" + c.c.](6.2) 

where the U{1) x U{1) symmetry (for A = 0) of the dual 
action has been elevated to a U{1)'^^ invariance (un- 
der independent phase rotations of each tpi field). Of 
this, only the single U(l) subgroup of identical rotations 
of all 2N fields has been gauged with the non-compact 
gauge field A^^. The A term breaks the C/(l)^^ symme- 
try to C/(l)^, of which iV — 1 are global, and the single 
gauge U(l) is preserved. In addition, there is a residual 
global symmetry under ^pj e^^-i'^/^tpj, ipj+N 
e-™j'^/Vi+Ar, with n, G {0, 1, 2, 3}, for j = 1 ... TV. 

In the N = 00 limit, the theory may be analyzed by 
saddle point methods. In particular, consider for sim- 
plicity the partition function with A = 0, which may be 
formally written 



Z 



where 
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— In / [dip] exp< — / d^rdr 



|(9^-a^)^|2+r|7/-|2 + i.|V'|4 \ (6.4) 



is the effective action of the D = i XY model as a func- 
tional of Afj_. Formally, at = 00, from Eq. 16.31 a 
saddle point approximation in A^ is justified, with the 
saddle point value being zero, A*^ = 0. The relevance 
of A is then determined by the two-point function of 

'^v\s:N = YH=i{i'*i't+N)'^ in the saddle point theory 
with A^ — 0. Since in this theory the ipi are decoupled 
XY fields (fiuctuating according to the non-trivial ?>D 
XY fixed point), one has then 
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A=0 



(6.5) 



N\{{r{x))\i^{x')f) 



N 



\x-x'\^^'^' 



The final expression obtains since the expectation value 
in the second line is none other but the two-point function 
of the four-fold symmetry breaking field at the D — i XY 
fixed point. Hence, one has A = 2A4, implying A > 6 
since A4 > 3. Thus in this limit the "monopole" (sym- 
metry breaking) terms are strongly irrelevant. We note 
in passing that the irrelevance of A can also be established 
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by working for A 7^ 0, and taking the large N hmit. The 
saddle point remains at = 0, and irrelevance follows 
simply from the irrelevance of "bi-quartic" coupling be- 
tween two 3D XY models by their four-fold symmetry 
breaking fields. 

We conclude that a direct second order transition with 
irrelevant instanton tunneling events is possible in this 
easy plane case. Note the crucial role played by the Berry 
phase term for the instantons in reaching this conclusion. 
Indeed it was the Berry phases that forced quadrupling 
of instantons thereby increasing their scaling dimension 
and making it possible for them to be irrelevant. 

While the A term may be irrelevant at the critical fixed 
point it is clearly very important in deciding the fate of 
either phase. In particular in the paramagnetic phase it 
picks out the particular pattern of translation symmetry 
breaking (columnar versus plaquette) and forces linear 
confinement of spinous. In critical phenomena parlance, 
it may be described as a dangerously irrelevant pertur- 
bation. 



B. Isotropic magnets 

In the context of the SJ models, the results of previ- 
ous sections show that the iV = 1, = 00, and easy 
plane N = 2 models all provide the same picture. A di- 
rect second order transition between the z condensed and 
VBS phases is possible with a 'deconfined' critical point. 
Right at this point, monopole tunneling events become 
irrelevant and spinon degrees of freedom emerge as the 
natural fields of the critical theory. This provides strong 
evidence that the same thing happens for the SU (2) sym- 
metric model {i.e at N = 2). 

What then is the proposed description of the critical 
point in the SU{2) symmetric model? This is simply the 
CP^ model with a non-compact gauge field and no Berry 
phase terms in Eq. (|1.7|) . Equivalently it may be thought 
of as the critical point of the D — 3 classical 0(3) model 
when monopoles have been forbidden by hand. This 
transition was first studied by Kamal and Murthy'^^ and 
more recently by Motrunich and Vishwanath'^^, where it 
was established that a continuous transition indeed ex- 
ists that is different from the Heisenberg transition. The 
non-compact CP^ theory Eqn. (|1.7II was also directly 
studied via numerical Monte Carlo methods and found 
to possess a continuous transition with the same univer- 
sal properties as the monopole suppressed 0(3) NLctM. 
Numerical results for exponents associated with several 
observables are available. Further evidence for the con- 
tinuous nature of the transition in the CP^ model coupled 
to a noncompact gauge field is obtained by considering 
the larger class of models with CP^~^ fields coupled to 
a noncompact gauge theory. It is well known that the 
N — 1 model has a continuous transitior^ which is dual 
to the XY transition, and a continuous transition is also 
expected for large values of N. Thus, the model of in- 
terest TV = 2 is sandwiched between these two extremes 



where a continuous transition is well known to obtain. 



VII. PHYSICAL PROPERTIES NEAR THE 
'DECONFINED' CRITICAL POINT 

We now discuss the consequences of the theory for the 
physical properties near the direct Neel-VBS transition. 
We will first discuss those properties that follow generally 
from the (dangerous) irrelevance of monopoles. Later 
we will specialize to the easy plane limit where the self- 
duality enables more progress. 

It is useful to think first about the various length scales 
in the problem in the VBS phase. 

First there is the spin-spin correlation length ^ which 
will diverge at the transition. Second, there is a length 
scale ^vBS associated with the 'thickness' of the domain 
walls of the (discrete) VBS order. The latter is clearly 
determined by the strength of the quadrupled monopole 
operator, A = A4, in Eq. H1.9|l : in the easy plane case, 
A appears as the co-efficient of a local term in the dual 
action Eq. H5.6|l . These two length scales will diverge 
differently - the domain wall thickness will diverge faster 
than the spin-spin correlation length. One can determine 
the scaling of ^vbs with ^ by a matching argument. On 
scaling grounds, one expects 



evBS - e/(AC'-^), 



(7.1) 



where / is a scaling function, and 3 — A is the RG eigen- 
value of A assuming the scaling dimension of the four 
monopole operator is A > 3, and d = 2, z = 1. Be- 
yond the scale of the correlation length ^, one can re- 
gard the VBS phase as XY ordered in 4'vbs, though with 
very weak four-fold anisotropy since A is irrelevant at the 
QCP. Hence the low-energy variations of the phase 9 of 
V'VBS ~ IV'VBsle'^ are described as a pseudo-Goldstone 
mode, with energy 



E{e) = / d^x 



y|V6i|2-Acos46' 



(7.2) 



where K and A oc A are renormalized parameters on the 
scale of ^. A twist of 9 (of e.g. tt/2) is carried hence by 
a domain wall which, by dimensional analysis, has width 

^VBS ^ \J K jX. Knowing then that ^vbs ^ A^^/^, one 
requires that f{x) ~ x~^/''^ in Eq. (|7.1|) . which implies^ 



VBS 



(A-l)/2 



(7.3) 



Since (A — l)/2 > 1, ^vbs indeed grows more rapidly 
than ^ as the QCP is approached. 

Thus there are two independent diverging length 
scales. Either of these length scales may be given several 
different interpretations. For instance, the spin correla- 
tion length ^ may also be interpreted as the length scale 
at which correlations of the dual global order parameter 
crossover from that of the critical fixed point to that of 
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the (unstable) fixed point which breaks the dual global 
continuous symmetry. Similarly, the domain wall thick- 
ness ^vBS of the VBS order is also the length scale at 
which the photon that couples to the spinous acquires a 
mass due to instanton effects. This is also the length be- 
yond which the logarithmic Coulomb potential between 
spinous crosses over to a linear confining one. This is 
distinct from the 'confinement' length scale describing 
the spatial size of the resulting two-spinon bound states 
(triplons) . This length scale is actually a non-trivial com- 
bination of the two other diverging length scales. It how- 
ever diverges faster than the spin correlation length 

Note that the critical theory is isotropic in space-time 
and therefore has dynamic scaling exponent z = 1. The 
values of other critical exponents may be obtained from 
the numerical work of Ref. 0. In the 0(3) symmetric 
case the correlation length exponent v ^ 1 (for and 
the Neel order parameter exponent (3 « 0.80. Perhaps 
most remarkably the anamolous dimension of the Neel 
order parameter field rj is large (« 0.6). This should 
be contrasted with the extremely small value for r/ at 
the usual Wilson-Fisher 0(3) fixed point in D = 3 di- 
mensions (and indeed for most other familiar three di- 
mensional critical points). The large value of rj can be 
rationalized by the thinking that the Neel order param- 
eter field decays into spinous right at the critical point. 
Indeed as argued in previous sections it is the spinous 
which appear as the more natural degrees of freedom at 
the deconfined critical point. We note however that the 
spinous are not to be considered 'free particles' - they are 
critical and furthermore interact through the coupling to 
the non-compact gauge field. 

Consider the effect of twisting the boundary condi- 
tions on the VBS order - for instance, for columnar 
dimerization prefer even columns at one boundary and 
odd columns at the opposite boundary. Let us suppose 
the twist is applied between the top and bottom ends 
{y = 0, W) of an L X W sample. On general grounds, 
the energy cost at long scales will he E ^ aL'^^^ = aL, 
where tr is a "surface tension" or domain wall energy 
per unit length. This surface tension is set, however, by 
the irrelevant monopole term and vanishes in a manner 
set by the divergence of the domain wall thickness. In 
particular, the surface tension scaling obtains only for 
twists sustained over a distance W > ^vbs- For twists 
of the VBS order over a shorter distance W such that 
^VBS ^ W ^ ^, the energy cost for this twist is greatly 
reduced to E ^ KL/{2W), where K is the "stiffness" as- 
sociated with the continuous dual global symmetry (and 
we are at length scales where the system has not realized 
this symmetry is actually discrete). The two energy costs 
for the twist become comparable for W ~ £,vbs, so that 
one expects a ~ K/S,vbs- This stiffness itself vanishes 
upon approaching the quantum critical point in a man- 
ner set by the divergence of the spin correlation length. 
Furthermore, the corresponding exponent is the same as 
for the spin stiffness on the other side of the transition. 
Specifically, the VBS stiffness K ~ ^'^-d-z ^-^^j-e ^ is 



the dual correlation length, d = 2 is the spatial dimen- 
sion, and z = 1 is the dynamic critical exponent. Thus 
K (X 1/C 

Note that this is not a test of self-duality but rather 
a test of the irrelevance of monopoles: the scaling of the 
VBS stiffness is a consequence of dual current conserva- 
tion which obtains if monopoles are irrelevant. Thus the 
same behavior is also expected for the isotropic model. 

In practice, a measurement of the domain wall energy 
in the columnar state is likely best obtained by simply 
comparing energies of systems of size LxW (in the x and 
y directions, respectively) with periodic boundary condi- 
tions in both directions {i.e. on the torus) and varying 
W. In particular, let us consider W > L, with L odd. In 
this case, the columns will prefer to align along the short 
direction (i.e. columns parallel to the x axis, breaking 
translational symmetry along y) in order to avoid intro- 
ducing a domain wall (which would be required in the 
other orientation due to the odd L) with energy cost aW. 
If W is odd, there will still be a domain wall required, 
but it is shorter and less costly, with energy aL. Hence 
one expects the ground state energy of the system to be 
Elw ~ eWL for W even, and E^w ~ eWL + aL for W 
odd. Here e is the ground state energy density. Hence 
the surface tension can be obtained by 

El.w+1 - ^Elm + El.w-1 ~ 2aL(-l)^, L odd. 

(7.4) 

This behavior will obtain provided W > L ■fvBS- 
When the system is smaller than the domain wall thick- 
ness, however, this energy is determined instead by the 
stiffness, i.e. 

Elm+1 - 2El,w + El^w-1 - K^i-l)"^, L odd, 

(7.5) 

for ^ < L < Vl/ < fvBS- 

It is also clear that at the critical point, both columnar 
dimer and plaquette order parameters will have power 
law correlators with the same exponent. This is indepen- 
dent of which one of these two phases we eventually end 
up in. This is because both order parameters are con- 
tained in the dual global boson creation operator. They 
correspond to the phase of this boson locking in different 
directions. Which one of these phases is selected is de- 
termined by the sign of the anisotropy term, but as this 
is renormalizing to zero at the critical fixed point, there 
will be no distinction between the two order parameter 
correlations at the critical point. 

This suggests the following interesting numerical 
check. Consider, for instance, the situation where the 
disordered phase has columnar VBS order. Now consider 
measuring the stiffness to plaquette order in this colum- 
nar phase. Since this order is not spontaneous in the VBS 
phase, measuring this stiffness cannot be accomplished as 
above by simply comparing systems with odd and even 
lengths. Instead, one should imagine introducing e.g. 
two rows separated by half W on which the magnetic 
couplings have been increased or decreased in a pattern 
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mimicking the strong bonds of the plaquette state. The 
ground state energy of this system should be compared 
with that obtained by shifting one of these rows by one 
lattice spacing, and the difference of these two energies 
interpreted as the energy cost for a twist of the plaque- 
tte order parameter. Deep in the phase, this energy cost 
will be exponentially small, AiS ~ Le"^/^^^^. However, 
in the regime ^vbs '> W ^ the cost for the twist of 
the plaquette order parameter will be determined by the 
stiffness, i.e. AE ~ KpiL/{2W). There wiU thus be a 
dramatic change from exponential to power-law behavior 
in this quantity on approaching the critical point. 

The coefficient of proportionality Kpi will equal the 
corresponding stiffness K for columnar order and scale 
identically to the physical spin stiffness on the other side 
of the transition. Once again, this is not a test of self- 
duality but follows from the irrelevance of monopoles and 
will hold in the isotropic case as well. 

Specializing to the easy plane case, the self-duality of 
the critical fixed point implies some further interesting 
properties. First, it is clear that the dual global boson 
will have the same power law decay as the physical spin- 
correlator. The former is identified with the VBS order 
parameter while the latter corresponds to the staggered 
XY correlators in the microscopic spin model. We thus 
have the remarkable result that the columnar dimer, pla- 
quette, and staggered XY magnetization all decay with 
the same power law right at the critical point. Further- 
more, the P exponent for the particular VBS order that 
actually develops will be the same as the /3 for the spin 
order. This is because the anisotropy only serves to lock 
the phase of the dual order parameter. The amplitude 
is already non-zero in the scaling limit near the critical 
point. 



A. Ordered state 

As elaborated in the previous subsection, there are two 
diverging length (or equivalently time) scales upon ap- 
proaching the transition from the VBS side. How does 
this manifest itself in the Neel ordered side? To under- 
stand this first note that in the ordered phase close to the 
transition there will be 'soft' modes that correspond to 
the incipient VBS order on the other side of the transi- 
tion. Indeed the frequency of these modes will go to zero 
at the critical point. For concreteness consider the case 
where the VBS order that develops is columnar. Then 
as is natural there will be a soft mode corresponding to 
columnar order with frequency vanishing on approaching 
the critical point. Remarkably there will be other soft 
modes that correspond to plaquette ordering whose fre- 
quency also vanish on approaching the transition. Once 
again this is despite there being no such order in the VBS 
phase. This result is already implied by the discussion 
in the preceding subsection. Indeed both the plaquette 
and columnar order parameters have power law correla- 
tions at the critical point regardless of which one of the 



two orders actually develops in the VBS phase. Thus 
it is natural that the frequencies of both modes go to 
zero on approaching the critical point from the ordered 
side. Furthermore, both soft mode frequencies vanish 
in exactly the same universal way on approaching the 
transition. Formally the columnar and plaquette order 
parameters are distinguished only by the orientation of 
the phase of the complex VBS order parameter. In the 
absence of monopole tunneling events (which tend to pin 
this phase to certain values), these two distinct order 
parameter fields will both behave in a common manner 
determined by the complex VBS order parameter. Thus 
in the scaling limit near the critical point both the pla- 
quette and columnar order parameters will display the 
same universal behavior. On general scaling grounds we 
expect the VBS soft mode frequency wvbs ~ Ps where 
is the ground state spin stiffness of the Neel ordered 
state. 

Despite the common universal behavior of the vanish- 
ing frequency of the two distinct VBS soft modes, there 
will be a small splitting between the two frequencies that 
is due to the irrelevant (quadrupled) monopole tunnel- 
ing events. Indeed the information about which of the 
two VBS orders eventually develops in the paramagnetic 
phase is contained in this small splitting. If columnar or- 
der develops, then the corresponding soft mode will have 
slightly lower frequency in the ordered phase as com- 
pared to the plaquette soft mode. Note that this split- 
ting will go to zero as the critical point is approached 
as the monopole fugacity renormalizes to zero. This will 
however vanish in a very different way from the overall 
VBS energy scale Wvbs- Indeed this splitting defines a 
new energy scale that vanishes faster than wyes- Thus 
we see that two different energy scales also characterize 
the physics of the ordered phase. 

In the easy plane case, one intriguing aspect of our the- 
ory is the physics of the vortex cores in the XY ordered 
phase close to the transition. As discussed extensively, 
there are two kinds of classical meron vortices which tun- 
nel into each other in the quantum theory. However, the 
irrelevance of these instanton tunneling events near the 
transition implies that the Ising order in the core will 
survive for a very long time. 

We have, up to this point, not considered any effects 
which explicitly break the lattice translational symme- 
try. Hence the discussion should be read as appropri- 
ate for extended, plane wave states of vortices. Cru- 
cially, these states may be classified by their (quasi- 
)wavevectors. Consider the continuum theory in which 
the spatially-oscillating single instanton events are ne- 
glected, and only the (irrelevant) quadrupled instanton 
fugacity A is included. In the dual formulation, the two 
Ising vortex states appear as relativistic charged parti- 
cles. They carry a conserved U(l) non-compact gauge 
charge (physically their vorticity), i.e. the number of 
these vorticies A''i -|- N2 is conserved. There is a dis- 
crete global symmetry, which implies that Ni — N2 
is conserved modulo 8. The latter is a consequence of 
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the continuum limit which removes the single instanton 
events, and is promoted to a continuous U(l) symmetry 
(with Ni — N2 fully conserved) if the quadrupled instan- 
tons are also neglected. Physically, then, excitations of 
the vortex vacuum, i.e. the XY ordered state, can be 
labelled by these quantum numbers. And one should ex- 
pect there to be "quasiparticles" (really "quasivortices" 
or "quasimerons" ) carrying an elementary unit (±1) non- 
compact gauge charge (physically ±2tt vorticity) and an 
elementary Z4 charge Ni — N2 = ±1. These vortices (as 
befits an XY ordered phase), with a 'core energy' scaling 
as 1/^. There will be other gapped excitations (also with 
a gap of 0(1/^)) with zero non-compact gauge charge 
(zero vorticity) and Ni — N2 = ±2, which can be viewed 
as the VBS soft modes, or alternatively as "vortex exci- 
tons" . 

Let us now consider the effects of the neglected oscil- 
lating single monopole term. Naively, this violates the Z4 
conservation law, and can mix the quasimeron states with 
Ni — N2 — ±1. However, due to the 4-sublattice oscilla- 
tion, the process in which one meron state is converted 
to another is accompanied by the addition of a large mo- 
mentum. In the absence of a sink for this momentum, 
therefore, even the single instanton term cannot mix the 
two meron quasiparticle states. At non-zero temperature, 
thermally excitated excitations with a large gap (3> 1 /C) 
can provide this momentum, but are present only with an 
exponentially small probability due to their gap. Hence 
violations of the conservation of the "Ising" charge of a 
single quasiparticle are exponentially weak at low tem- 
perature. Of course, four quasiparticles with "up" Ising 
cores can scatter off one another to produce four quasi- 
particles with "down" Ising cores via the non-oscillating 
A term. The amplitude for even this process is, however, 
suppressed by the irrelevance of A if one is near the QCP. 

Thus we arrive at the remarkable conclusion that the 
"elementary" gapped vortex excitations of the XY or- 
dered phase carry a sharp extra Ising quantum number. 
Ramifications will be explored in Sec. IIX A II 

A number of other predictions may also be made on 
the effects of various perturbations that can be added to 
the Hamiltonian near the zero temperature critical point. 



B. Uniform Zeeman field 

First let us think about the excitation structure in the 
paramagnetic side. Deep in the phase, the lowest spin- 
ful excitations will be S" = 1 which will be gapful. On 
approaching the transition, due to the diverging 'confine- 
ment' length one might naively think that this will break 
up into spinous. However even with a non-compact gauge 
field due to the log attraction coming from the photon 
there will be logarithmic confinement and the lowest en- 
ergy spin carrying excitations will continue to have S = 1 
This "magnon" is a gauge neutral bound state of two 
spinous. Now imagine sitting in the paramagnetic side 
close to the transition, and turn on a Zeeman field along 



the z— axis in spin space. Once the Zeeman energy ex- 
ceeds the magnon gap, the chemical potential for such 
magnons becomes positive and they should condense to 
modify the ground state. This leads to XY antiferromag- 
nctic order. In the non-compact approximation, because 
this condensate is gauge neutral, it does not create a 
gap for the photon via the Higgs mechanism. In real- 
ity (i.e. beyond the non-compact approximation) what 
this means is that there will be coexistence between VBS 
and XY order. As the spinous are not condensed in this 
phase, there is no disruption of the VBS order (this can 
be seen e.g. from the fact that the spinous appear as 
dual "vortices" in the VBS order parameter, while the 
magnons carry zero dual vorticity). 

These considerations hold only provided another tran- 
sition does not pre-empt magnon condensation as the 
Zeeman field is increased. This will happen, e.g. if the 
magnons experience attractive interactions with one an- 
other. Indeed, in a Coulomb interacting system, it is 
natural to expect that the magnons, which are the ana- 
log of excitons, will have attractive interactions with one 
another at long distances, due to the analog of "van der 
Waals" forces between their fluctuating dipole moments. 
This attraction, however, competes with other local in- 
teractions due to the complex critical physics on scales 
^ ^, so the outcome is not clear. Therefore, we do not see 
a clear argument against a continuous magnon conden- 
sation transition into a coexistence phase. Likewise, of 
course, we cannot rule out a direct first order transition. 
In any case it therefore seems as though a direct second 
order transition between the magnetically ordered and 
VBS phases is unlikely at non-zero Zeeman field. 

If the coexistence phase exists, it is interesting to con- 
template the transition between the coexistence phase 
and the pure magnetically ordered one (with canted an- 
tiferromagnetic order). In the system with XY symme- 
try, a transition with exactly these symmetries has been 
studied in Ref. l56t In this paper, it was shown that, 
despite coexisting superfluid order and the consequent 
gapless goldstone mode, this transition is in the univer- 
sality class of a D = 3 XY model, the Z4 symmetry 
breaking perturbation and the coupling to the goldstone 
mode being irrelevant. Thinking in terms of the dual 
formulation suggests this analysis should apply here. In 
particular, both vortex fields ipi^2 remain gapped across 
the transition, since both phases are XY ordered. Only 
the composite order parameter i/jvbs = ^i'02 involved 
in the criticality, and it does not couple minimally to the 
non-compact gauge field. Integrating out the massive 
vortex fields while keeping a composite ^ybs field and 
the gauge fluctuations describing the Goldstone mode, 
one arrives at a model equivalent to Ref. Is^f Note that 
this result implies that this critical point is also decon- 
fined in precisely the same sense as the others discussed 
in this paper. 
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C. Staggered Zeeman field 

Consider the effect of a staggered Zeeman field on the 
original spin model. First assume easy plane anisotropy 
in the plane orthogonal to the applied field. The stag- 
gered field will always induce some staggered magnetiza- 
tion but we can ask about XY or VBS order superim- 
posed on this. In the CP^ description, a staggered Zee- 
man field corresponds to a uniform 'magnetic' field that 
couples to the z component of h. It is quite clear that 
there will now be a split transition between the VBS and 
XY ordered phases with an intermediate phase with nei- 
ther XY nor VBS order (but of course with a staggered 
magnetization) . 

Consider the approach from the VBS phase. In the 
presence of the staggered field, one of the two CP^ fields 
will condense first. This transition is described by the 
iV = 1 SJ model (and is inverted D = 3 XY). The re- 
sulting phase is the advertised phase with neither XY or 
SP order. Actually it is more useful to think of this crit- 
ical point in the gauge language as a 'deconfined' critical 
point than just as inverted XY. This is so particularly if 
one asks about the magnon spectral function at this tran- 
sition. This will be determined by the spinon dynamics 
which in turn are coupled to a non-compact U(l) gauge 
field. Thus we might again expect anamolously broad 
spectral functions even though both phases are confined. 
(Note that total is still conserved). 

Eventually, as one tunes towards the XY ordered 
phase, the other CP^ field will also condense leading to 
XY order. This transition is actually exactly dual to the 
other one discussed above. This is because the staggered 
Zeeman field couples to the same operator in both the 
CP^ and dual representations. For small staggered Zee- 
man fields, the phase boundaries must come in with the 
same exponents, etc. 

Finally, in the 0(3) model, the second transition will 
not happen as the Neel vector will immediately line up 
with the staggered field. However, the first transition will 
continue to be described by the TV = 1 SJ model. Details 
of the slopes of various phase boundaries etc. may be 
found Ref. |^ - in the terminology of that reference, a 
uniform field corresponds to the staggered Zeeman field 
discussed here. 



D. Finite temperature transitions 

Finite temperature properties near the transition may 
also be discussed. Here it is clearly necessary to distin- 
guish between the easy plane and isotropic cases. In the 
latter, the Neel order does not survive for any non-zero 
T while in the former case there is power law order at 
low T which eventually disappears through a Berezinski- 
Kosterlitz-Thouless (BKT) transition. In both cases 
however the discrete broken lattice symmetry of the VBS 
phase will survive upto a non-zero finite T . The associ- 
ated finite-T transition will be in the universality class of 



the .^4 clock model in cZ = 2. This transition is known to 
be described by a line of fixed points with continuously 
variable [i.e. nonuniversal) critical exponents. The line 
of fixed points results from the exact marginality of the 
fourfold symmetry breaking term V'vbS' ^^"^ coeffi- 
cient. As A approaches zero, the nature of the criticality 
approaches that of a simple D = 2 classical XY model, 
i.e it becomes BKT-like. Thus i] ^ 1/4 and f ^ I/|A| 
diverges in this limit. Since at the zero temperature QCP 
instantons are irrelevant, we may conclude that the fixed 
point value of A at the classical VBS-paramagnet tran- 
sition (which is generally finite and non-zero at Tc) ap- 
proaches zero as Tc —> 0. Hence the non-universal critical 
behavior of the VBS-paramagnet transition becomes ar- 
bitrarily close to BKT behavior as this transition line is 
followed into the T — QCP. This conclusion is inde- 
pendent of the XY or 0(3) symmetry of the magnetic 
ordering. 

In the XY case, the self-duality of the easy plane fixed 
point implies further that the phase boundaries associ- 
ated with the finite-T transition from both the Neel and 
VBS phases have the same shape at low T. Note that 
both are BKT-like for Tc asymptotically close to zero 
{i.e. near the QCP), consistent with duality. Indeed, one 
expects not only the phase boundary but also all critical 
correlations to match in this limit, including amplitude 
ratios. 



VIII. DECONFINED QUANTUM CRITICALITY 
AT THE VBS TO SPIN LIQUID TRANSITION 

In this section we argue that the transition between a 
valence bond solid and a fractionalized spin liquid is also 
an example of a deconfined quantum critical point in a 
precise sense. 

In two spatial dimensions, a fractionalized spin liquid 
is expected to be described as the deconfined phase of 
a Z2 gauge theory with a gapped Z2 vortex - the vi- 
son. This Z2 gauge field is minimally coupled to spin-1/2 
spinon excitations. We only consider the case where the 
spinous are gapped. A precise theoretical characteriza- 
tion is given by the notion of topological orderiiSiS. 

Consider the evolution of the ground state of a spin- 
1/2 system (or equivalently for bosons at half-filling) be- 
tween such a fractionalized spin liqiud and a VBS on, say, 
a square lattice. Despite the lack of any obvious local or- 
der parameter for the spin liquid there is a close similar- 
ity with the Neel- VBS transition. Indeed both the spin 
liquid and VBS are characterized by two distinct types 
of order (the former by topological order and the latter 
by broken lattice symmetry) . Naive thinking might then 
suggest that a direct second order transition is not possi- 
ble. Rather one might have expected two transitions with 
an intermediate "coexistence region" which breaks lat- 
tice symmetry but is also topologically ordered (a VBS* 
phase, in the notation of Ref. ^3)- Once again this naive 
expectation is incorrect and a direct second order tran- 
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sition is indeed possible. Furthermore the critical theory 
may be regarded as a non-compact U(l) gauge theory 
with an extra emergent dual global U(l) symmetry. 

It is convenient to begin with a theoretical formula- 
tion that is powerful enough to describe both phases and 
all of their distinct excitations. Such a formulation is 
provided in the work of Refs. 'l5"l()"58'. As before, the 
underlying spin model is first reformulated as a theory 
of spin-1/2 spinon fields that are minimally coupled to a 
compact U(l) gauge field with Berry phases. The VBS 
corresponds to a confined paramagnet where the spinous 
have disappeared from the spectrum. The spin liquid 
obtains when a singlet pair of spinous - which carries 
gauge charge-2 - condenses i.e enters a Higgs phase. Let 
us represent this Higgs field by the operator Q — e*'''. 
We imagine integrating out the individual spinon fields 
(this is permissible because all spin carrying excitations 
are gapped in both the VBS and spin liquid phases) , and 
obtain the following theory for the transitionpii^ 

S = Sip + Sa + Sb 

= -2i^cos(A</3-2a), (8.1) 

e 

where Sa and Sb are defined in Eqs. H3.12|l and H3.13|l . 
Note that Eq. (|8.1|) is just the = 1 SJ model studied in 
Section FlV Al but with the crucial difference that ip car- 
ries charge 2 (compare with Eq. H4.1|l 'l. The duality trans- 
formations of Section llV Al are easily applied to Eq. H8.1|l . 
and we obtain an XY model with 8-fold anisotropy which 
is irrelevant at the transition (this contrasts with the 4- 
fold anisotropy obtained in Section FlVAp . 

In the spirit of previous Sections (particularly Sec- 
tion ^^3, these results may be understood physically 
as follows. In the fractionalized phase the condensation 
of the charge-2 scalar leads to vortex excitations (the 
visons) which carry tt gauge flux. In the fractionalized 
phase instanton effects kill visons in pairs - indeed this is 
precisely what leads to their 'Ising' nature. The transi- 
tion to the confined VBS phase occurs when the visons 
condense. But near the transition, and in the continuum 
limit, we expect once again that all instanton events are 
quadrupled. Thus the tt flux vortices can only disap- 
pear 8 at a time. This gives the XY model with 8-fold 
anisotropy. 

We thus see that the 8-fold anisotropy in the dual 
XY model should be interpreted as instanton tunnel- 
ing events in the original compact gauge theory. Con- 
sequently as before we conclude that instantons are ir- 
relevant at the critical fixed point so that a gapless non- 
compact U(l) gauge theory obtains. (We remind the 
reader that the global XY model is the dual of the con- 
densing charge-2 scalar coupled to a non-compact gauge 
field). 

Note once again the crucial role played by the Berry 
phases which are responsible for leading to an 8-fold 
anisotropy (as opposed to 2-fold as would obtain in their 
absence). 



Note also that spinous are well-defined in the fraction- 
alized phase but are confined in the VBS phase. What 
is the fate of the gapped spinous right at the transition 
between these two phases? The arguments above show 
that at the critical point the spinous are minimally cou- 
pled to a non-compact U(l) gauge field descending from 
flp (which in turn is also coupled to the critical spinon 
pair field) . The strong scaling properties of this critical 
point (which is dual to the D = 3 XY model) implies that 
the gauge-field has the following two-point correlator at 
criticality 

(a,ip)aA-p))--(s,.-^), (8.2) 
p \ p^ J 

where is the spacetime 3-momentum. Note that this 
propagator does not have the Maxwell scaling, but a 
l/p dependence fixed by the scaling dimension dim[a^] — 
1. This implies a l/r interaction between static massive 
spinous at criticality. 

It is sometimes stated that the transition between the 
VBS and spin liquid phases is described by a Z2 gauge 
theory. The results here however show that the transition 
is in fact described as a deconfined U(l) gauge theory in 
a very precise sense. It is the spin liquid phase itself (as 
opposed to the transition) that is described as a (decon- 
fined) Z2 gauge theory. 

A. Spin liquids that break lattice symmetry 

An important subtlety has been glossed over in the 
analysis so far in this section. Spin liquid states with no 
broken lattice symmetries are certainly possibleiS*i2i^2iSfij 
and for these the above analysis applies. However, in 
bosonic mean field theories of SU(2) spin liquid states 
on a variety of lattices^^iSiiSSiS^iSij the spin liquid state 
is commonly found to break a global lattice rotation 
symmetry^^ - such a state has 'bond-nematic' order. The 
spin liquid is associated with short-range, incommensu- 
rate spin correlations at a wavevector K , and the choice 
of K often breaks a lattice symmetry e.g. a spin liq- 
uid state at K = (fc, k) is distinct and inequivalent to a 
state at K — {k,—k). Such states appear naturally at 
the boundary of a VBS stated, and for these the the- 
ory above has to be reconsidered. Before doing this, 
we note one important case for which this is not nec- 
essary: the CS2CUCI4 latticeSLS, which interpolates be- 
tween the square and triangular lattices. Within a large 
A^ bosonic mean field theory treatment, the ground state 
in the square lattice limit is a VBS, while in the triangular 
limit it is spin liquid which breaks no lattice symmetries: 
the transition between these states is described by the 
theory in Eq. H8.1|l . 

Turning to a spin liquid that does break lattice symme- 
tries, consider e.g. the transition on the square lattice^ 
from the VBS (Fig ^| with short-ranged spin correla- 
tions peaked at the wavevector (tt, tt), to a bond-nematic 
spin liquid at wavevector K = {k,k) or K = {k,—k). 
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The choice of either of the latter states breaks a symme- 
try of reflection about the principal square lattice axes. 
In mean field theoryj^ this transition is characterized 
by the condensation of two Higgs fields, which we de- 
note as Qx — e"^^ and Qy = e"^". These fields are 
odd under the lattice reflections TZx and TZy in Table HI 
respectively,^^, and this prohibits terms which are linear 
in either Higgs field in the effective action. Using these 
symmetries, and the requirements of gauge invariance, 
we generalize Eq. (I8.I|I to 

S = Sip + Sa+ Sb 



Sep = — |^2i cos {A.ipx — 2a) 

e 

+2t cos {A(py — 2a) + 2t' cos(2<y92, — 2(py 



(8.3) 



where, as before, Sa and Sb are defined in Eqs. H3.I2|I and 
(|3.I3I) . Note the crucial factor of 2 in the argument of the 
third cosine in Eq. (|8.3(l : this is required by the inversion 
constraints above. Apart from the usual compact U(l) 
gauge invariance, Eq. (|8.3|l is also invariant under the 
global Z2 transformation 



fx 'Px+ 7r/2 
ipy ~* ipy^ n/2 



(8.4) 



which realizes the lattice reflection symmetry. (Note that 
the square of the transformation in Eq. (|8.4|) is equivalent 
to the identity modulo a compact U(I) gauge transfor- 
mation.) Consequently, there are now two inequivalent 
Higgs phases, with {tp^ - (fy) = (or (Q^) = (Qy)) and 
{ipx~ipy) — TT (or (Qx) — ~{Qy)), and these correspond^ 
to the two possible spin liquid phases at K = {k,k) and 
K = {k,—k). The theory Eq. H8.3(l can be analyzed by 
the same duality transformation applied to Eq. 1)8. I|l . but 
the critical properties have not been determined. 



IX. ANALOGIES AND EXTENSIONS 

A. Superfluid-insulator transition of correlated 
bosons 

The models and the phenomena discussed in this paper 
can be fruitfully discussed from a different point of view. 
Consider a system of bosons with short-ranged repulsive 
interactions on a square lattice such that there is half a 
boson per site on average. It has long been appreciated 
that such a bosonic system is closely related to quantum 
spin models with easy plane (or easy axis) anisotropy. 
Indeed, there is an exact equivalence in the hard-core 
limit in which at most one boson occupies each lattice 
site. Specifically, one may consider a model of bosons 
(described as 0(2) quantum rotors) on a square lattice: 



-i^cos( 



Here (pr G [0, 2tt) represents the boson phase, rir is the 
conjugate boson number and is an integer S [— oo,cx3]. 
The ellipses represent other short ranged terms that can 
be tuned to drive transitions from a superfluid to (for 
instance) the bond stripe insulator. To relate the above 
boson Hamiltonian to the antiferromagnetic systems con- 
sidered in the bulk of the paper, we note that one may 
define 



= ere 

S„ - — - 



T^4 



(9.2) 
(9.3) 



</>;) + .... (9.1) 



For large U this gives a faithful representation of an easy- 
plane spin- 1/2 antiferromagnet, and the universal physics 
is expected to be unchanged at smaller U. 

Clearly a superfluid phase of the bosons is possible 
(and corresponds to the XY ordered phase in the magnet 
analogy) . Various kinds of Mott insulating ground states 
are also possible. (These correspond to quantum param- 
agnets in the magnetic case). A simple Mott state cor- 
responds to the bosons forming a checkerboard ordered 
pattern in which the sites of one sublattice are preferen- 
tially occupied. This will be stabilized by large nearest 
neighbor repulsion and corresponds to the Ising ordered 
antiferromagnet. In the boson language the columnar 
VBS state may be understood as a 'bond-centered' stripe 
(or a bond density wave) - a state in which each boson 
is shared in a bond between two nearest neighbor sites 
such that these favored bonds have lined up in columns. 
The considerations of Section more or less apply di- 
rectly to the transition between the superfluid and the 
bond-centered stripe insulator (or the analogous plaque- 
tte ordered insulator). In particular the critical theory is 
'deconfined' and is expressed in terms of two fields each 
with boson charge-1/2 that are minimally coupled to a 
non-compact U(l) gauge field. However the discussion in 
Section |V| was intended for weak easy plane anisotropy 
on an isotropic spin model. It is somewhat more satisfy- 
ing to derive the crucial field theory Eq. H5.6|l directly for 
the bosonic system. We point out that the approach of 
Ref. provides such a direct derivation of the required 
dual action. However the close connection with fraction- 
alized charge degrees of freedom is somewhat obscured 
by that approach. We therefore sketch in Appendix |E| 
a derivation proceeding in a manner more similar to the 
considerations of the previous sections, in particular go- 
ing directly from the boson model of Eq. H9.1|l to the dual 
meron action obtained earlier. 

In the context of boson models (in view of potential 
applications e.g. to atomic bosons in optical lattices or 
to electronic systems where the bosons are Cooper pairs), 
some physical properties arise which are less natural in 
the context of quantum antiferromagnets discussed ear- 
lier. In particular, it is interesting to consider the effects 
of an applied orbital magnetic field coupling. This can 
bring out the unusual physics of Ising ordering in the 
vortex cores discussed earlier in Sec. IVII Al 
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1. Orbital magnetic field 

Let us consider the structure in an applied orbital mag- 
netic field B. The QCP at zero field describes a transition 
between a superfluid phase and a bond-centered striped 
phase. 

Suppose the system is on the superfluid (XY ordered) 
side of the QCP, and a small magnetic field is applied (we 
use the internal field B). This field produces vortices, 
separated by an average distance £ — yJcjiQ/ B, where 
4>Q = hc/q is the flux quantum, and q is the boson charge. 
In the weak field limit, where the length I is large, one 
expects these vortices to form an Abrikosov lattice, since 
the long-range logarithmic interactions between vortices 
dominate their kinetic energy. Now suppose one is near 
the QCP, so that the correlation length ^ is large. To 
a first approximation, one can neglect instanton events, 
and treat the Ising quantum number of the vortices as 
conserved. Then each vortex in the Abrikosov lattice has 
a definite Ising "charge" , and hence the system as a whole 
some sort of Ising magnetic order. It is straightforward to 
see that the basic interactions between these Ising "spins" 
are antiferromagnetic, and that these interactions decay 
rapidly when the two vortices in question are separated 
by a distance much larger than ^. These interactions 
arise because the two types of vortices carry opposite 
gauge fiux Jd?x eijdiUj = ±7r. This gauge flux is con- 
flned to a region of the size of the (gauge) "penetration 
depth" . Since the gauge field fluctuations are part of the 
quantum critical theory, this penetration depth is of the 
0{£,). Due to the Maxwell term in the action, two nearby 
vortices have lowest energy with opposite gauge fluxes 
(and hence smaller total gauge flux), provided the two 
fluxes overlap. A mean fleld analysis following Abrikosov 
leads to the same conclusion as discussed below. In par- 
ticular, consider the Lagrangian C{za) (in Eon. n~7|) for 
the Za fields. We are interested in w < (in Ean. I1.8|l . 
and it is convenient to consider the limit w = —2u + Sw, 
with < (5w <C 2u. For Sw — 0, the mean field the- 
ory (which neglects fiuctuations of a) comprises simply 
of two decoupled copies of Abrikosov's lattices for zi and 
Z2- Thus the solution consists of a triangular vortex lat- 
tice in each Za, with lattice spacing y/2£ (since each Za 
has charge q/2). These two lattices are completely de- 
coupled in this approximation. With Sw > 0, the energy 
is minimized when the integral of |2:ip|2;2p is smallest. 
This is accomplished by placing the zi and Z2 vortices 
as far apart as possible, so that |zip|z2p is reduced over 
the maximum spatial area. The solution is to choose the 
two triangular Za vortex lattices as the two sublattices 
of a honeycomb lattice. This corresponds to an antifer- 
romagnetic orientation of the Ising vortex cores on the 
honeycomb. 

To establish the stability of this order, we must recon- 
sider the effect of instanton events in this phase. The 
important events are single instantons, which act like a 
transverse field on the Ising quantum number. While 
these average away in the continuum theory, the finite 



lattice spacing £ provides an upper length cutoff for 
the oscillations of the single instanton fugacity, which can 
therefore have an effect. Near the QCP, it is possible for 
both £ and to be large, but to have £ not much greater 
than ^. In this limit, which we consider, the overlap of 
the vortex cores is strong, hence the Ising antiferromag- 
netic "exchange energy" between neighboring vortices is 
large, i.e. of order 1/^ by scaling. The effective trans- 
verse field on the vortices is more difficult to estimate. 
In a mean fleld treatment, one simply averages the oscil- 
lating instanton fugacity over the two dimensional ~ 
using some smooth envelope function. This gives a trans- 
verse fleld ^ Ao/^^. Fluctuation effects may be expected 
to further decrease this field. Hence, the transverse field 
is much weaker for large ^ than the antiferromagnetic 
coupling between cores. 

Thus we arrive at the remarkable conclusion that the 
vortex state near the QCP exhibits antiferromagnetic 
Ising LRO of the staggered Ising magnetization of the 
vortex cores. Note that this analysis applies when the 
magnetic length £ ^ For very small fields, or further 
from the QCP, £ ^ and the antiferromagnetic interac- 
tions between cores e~^/^) decay exponentially, while 
the transverse field is likely of power law form. Hence for 
very small fields it seems probable that the Ising cores 
become disordered. In this case, the physical manifes- 
tation of the long-lived Ising staggered magnetization is 
the presence of a low energy "antibonding" excitation of 
each vortex. 

Clearly, upon increasing the quantum fluctuations, this 
vortex lattice must disappear, since the VBS state on the 
other side of the critical point is an orbital paramagnet, 
i.e. no change in symmetry occurs on applying a weak 
fleld to it. Hence there is at least one phase boundary 
separating the Ising ordered vortex lattice from the VBS 
phase that persists at B > 0. We will not address this 
"vortex lattice melting" physics here, except to say that 
flrst order, continuous, and multi-stage transitions (with 
intermediate partially ordered phases) are all possible in 
principle (and difflcult to distinguish between on purely 
theoretical grounds). 



B. Higher spin 

In the bulk of this paper we have focused on the 
spin-1/2 square lattice antiferromagnetic model. Here 
we briefly discuss the fate of higher spin models on the 
square lattice. It should be clear by now that the answers 
depend crucially on the Haldane phases that obtain for 
higher spin. In the isotropic model, if 25* = l(mod 4) 
then the monopoles are quadrupled. Thus for all such 
values of the spin, a direct second order Neel-VBS tran- 
sition described by the same deconflned critical theory 
as for spin-1/2 obtains. If 25* = O(mod 4), then there 
are no oscillating phase factors for the monopoles. This 
has the consequence that a translation symmetric quan- 
tum paramagnetic state is now possible. The transition 



27 



to this state from the Neel state will be described by 
the usual LGW 0(3) fixed point {i.e. with monopoles 
present). If 25' = 2 (mod 4) then the appropriate Haldane 
phases lead to doubling of monopole events. Now con- 
fined paramagnetic states necessarily break lattice sym- 
metries. Whether a direct second order Neel-VBS tran- 
sition is allowed or not depends on the scaling dimension 
of the 2-monopole operator at the monopole-suppressed 
fixed point. If this is irrelevant, then the same deconfined 
critical theory as for spin- 1/2 will obtain. 

It is interesting to consider the spin-1 case in the pres- 
ence of some easy plane anisotropy. This may equiva- 
lently be viewed as a model of bosons at integer filling - 
unlike in the isotropic limit, a translation symmetric con- 
fined paramagnet is clearly possible. A direct transition 
between the XY ordered phase and such a paramagnet 
is clearly possible and will be in the usual D — 3 XY 
universality class. 

However, presumably the interesting question even in 
the easy plane case is whether a direct second order tran- 
sition is possible between the XY ordered phase and a lat- 
tice symmetry broken confined paramagnet with bond or- 
der. To answer this question and to obtain a description 
of such a paramagnet, it is convenient to start from the 
isotropic limit and introduce weak easy plane anisotropy. 
In the isotropic limit the confined paramagnetic states 
will break lattice symmetries and this will be preserved 
upon turning on some easy plane anisotropy. In a CP^ 
description, there will now be monopole Berry phases 
that oscillate on two sublattices of the dual lattice. We 
may now dualize in the easy plane limit to the meron 
vortex degrees of freedom. In this description it is clear 
that the translation broken VBS state is again described 
by an equal amplitude condensate of both vortex fields. 
Now the instanton term converts two merons of one kind 
into two of the other kind i.e. the co-efficient of the 
A term in Eq. is Re[(i/'tV'2)^], with A = A2. The 

relevance/irrelevance of this at the self-dual, easy plane, 
non-compact CP^ fixed point will determine whether a 
direct second order transition obtains or not: this ques- 
tion remains open at present. 

Note that near the transition to the usual paramag- 
net with no broken symmetries there will only be single 
species of vortex with a featureless 'paramagnetic' core. 
On the other hand near the transition to the VBS phase 
there will once again be two species of (nearly) stable 
vortices with cores that have very long-lived Ising order 
(if as discussed above such a direct continuous transition 
is possible). 



C. Honeycomb lattice 

The considerations of this paper generalize readily to 
other two dimensional bipartite lattices. For instance on 
the honeycomb lattice the Haldane Berry phase calcula- 
tion implies that all monopole events are tripled (rather 
than quadrupled as on the square lattice). This implies 



that the issue of whether or not a direct second order 
transition described by a 'deconfined' critical point ob- 
tains between the Neel and VBS phases is determined by 
the scaling dimension of the 3-monopole operator at the 
monopole-suppressed fixed point. Unlike the square lat- 
tice it is however less clear that the 3-monopole operator 
will be irrelevant. For instance for the TV = 1 SJ model, 
the Higgs-VBS transition is determined by the Z3 clock 
universality class which is distinct from the XY univer- 
sality class. In other words the 3-fold anisotropy which 
represents instantons is relevant at the deconfined fixed 
point at = 1. At large- A^ all monopoles continue to 
be irrelevant. The fate of the physical models (with and 
without easy plane anisotropy) can only be settled by di- 
rect numerical computation of the scaling dimension of 
the 3-monopole operator. 



D. Ising Anisotropy and Other Transitions 

In this paper, we have focused on the properties of 
spin-1/2 antiferromagnets with full SU(2) spin rotational 
symmetry, or its easy-plane reduction to U(l). These two 
cases are amenable to analysis due, on the one hand, to 
natural continuations of the SU(2) invariant CP^ rep- 
resentation to CP^^^, and through standard XY dual- 
ity. One may also ask whether similar deconfined critical 
points might arise in systems with easy-axis (i.e. Ising) 
anisotropy, which also retain the U(l) subgroup of SU(2). 
Unfortunately, this limit is much less amenable to micro- 
scopic duality transformations on the lattice level, and 
so it is difficult to make firm statements. While some of 
us suspect that no deconfined critical behavior is likely 
in this case, it is nevertheless of interest to present can- 
didate field theories for such deconfined transitions. 

Very naively, one may attempt to begin with the CP^ 
representation of the quantum antifcrromagnet, and sim- 
ply change the sign of the anisotropy term, taking w > 
in Eq. (|1.8|) . In a mean-field analysis of the continuum 
field theory of Eq. (|1.7|) , including anisotropy of this sign 
would indeed have the desired effect of yielding a tran- 
sition between an Ising ordered phase (for s < 0) and a 
VBS phase (for s > - actually one naively obtains the 
Coulomb phase of the gauge theory, neglecting the dan- 
gerously irrelevant instantons). There are, however, sev- 
eral caveats to this candidate theory that must be men- 
tioned. First, supposing the gauge field non-compact, 
fluctuation effects are known in some situations {e.g. the 
classical Abrikosov transition between the normal state 
and vortex lattice at Hc2) to drive naively continuous 
transitions involving gauge fields first order. While we 
believe this does not occur in the cases of SU(2) in- 
variance and XY anisotropy, these conclusions are based 
on s everal exact dualities and the numerical results of 
Ref. I21I39I directly on models in which instantons have 
been suppressed, and direct simulations of CP^ models 
coupled to non compact gauge fields-i. Second, even if 
the non-compact transition is continuous, to constitute a 
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stable deconfined QCP, it must be stable to the (quadru- 
pled) instanton events allowed by the microscopic com- 
pact model. 

At present we do not have supporting evidence in fa- 
vor of either of these two conditions. It would be of 
some interest to develop a semiclassical description of 
the above scenario to better evaluate it in physically in- 
tuitive terms. We note that XY anisotropy, for very sim- 
ple reasons, favors a "deconfined" critical scenario. In 
particular, weak XY anisotropy converts the topological 
defects (solitons) of the antiferromagnet from skyrmions 
to merons, "fractionalizing" them already in the antifer- 
romagnetic phase. Adding Ising anisotropy instead ren- 
ders the topological defects local "droplets" , or domains 
of antiphase ordering. These can apparently be viewed as 
distorted skyrmions, in which the smooth rotation from 
the anti-aligned core to infinity is replaced by a domain 
wall of finite width. Thus there is no fractionalization 
of the topological defects in the Ising antiferromagnet, 
although they do appear to carry the integer skyrmion 
quantum number. 

Nevertheless, the action in Eqs. H1.7|l and H1.8|l appears 
to describe a putative deconfined Ising AF to VBS tran- 
sition. There is clearly no self-dual description of this 
QCP, since neglecting instantons, the VBS phase is re- 
placed by a Coulomb phase with a gapless photon, while 
the Ising AF has no gapless excitations. Formally, how- 
ever, one may wonder what physics might be represented 
by considering the mathematically similar "anisotropy" 
in the dual meron theory, i.e. taking > in Eq. H5.6|l . 
Provided this transition remains continuous in 2 -I- 1 di- 
mensions and the A term remains irrelevant in this case, 
this would describe a different quantum phase transition. 
In particular, for Sd > 0, the ground state has no vor- 
tices and there is a Meissner response (Maxwell term for 
A^)^ hence it describes an XY superfluid. For Sd < 0, 
with w > 0, one or the other (not both) types of merons 
condense and the gauge field develops a Higgs mass. 
Hence this describes a non-superfluid state. From Ta- 
ble one sees that the non-zero expectation value of 
IV'iP ~ IV'aP implies Ising AF order. Thus this critical 
point describes a putative direct continuous transition 
between XY and Ising antiferromagnets. 

Thus these two theories describe different potential 
routes of "disordering" a quantum Ising antiferromag- 
net with U(l) spin-rotation symmetry, either to a VBS 
phase or an XY antiferromagnetic phase. If they are in- 
deed continuous, with the instanton fugacity and A term 
respectively irrelevant, they are not self-dual but instead 
dual to one another. It would be interesting to deter- 
mine with more certainty whether these putative critical 
theories can survive fluctuation effects. We note that 
the numerical simulations of Ref. |^ observed a first or- 
der transition between VBS and Ising antiferromagnetic 
phases. While this does not rule out the possibility of a 
continuous transition in other microscopic models, it is 
perhaps some evidence to the contrary. 



X. EXPERIMENTS 

We now briefly discuss the implications of the phe- 
nomena discussed in this paper for experiments on quan- 
tum magnets. What are some good signatures of these 
phenomena in experiments? Imagine a quasi two dimen- 
sional Mott insulator where each layer has a square lattice 
of localized spin-1/2 moments. Ignoring all effects due to 
coupling between the layers and to other degrees of free- 
dom (phonons,etc), a direct second order zero tempera- 
ture transition between the Neel ordered and translation 
broken VBS phases should be possible (for instance by 
application of pressure) with properties described as in 
previous sections. It is first important to emphasize that 
the proposal of a deconfined critical point is on firmest 
ground for a system with a spin-1/2 moment per unit cell. 
With higher spin or with more than one spin-1/2 moment 
per unit cell other (more ordinary) kinds of phase tran- 
sitions may well obtain. In practice (even with spin-1/2 
per unit cell) the growing VBS fluctuations associated 
with the lattice symmetry breaking will couple strongly 
to lattice disortions particularly at low temperature . If 
the phonons can be regarded as three dimensional (even 
though the magnetic interactions may be well approxi- 
mated as two dimensional), a small region of coexistence 
will most likely be introduced at very low temperature. 
This may be roughly understood as follows. The elastic 
energy cost of a latttice disortion of magnitude x that 
couples to the VBS order parameter is of order x'^ . How- 
ever the electronic energy gain is much bigger (as the 
susceptibility associated with the VBS order parameter 
diverges at the transition), going as x'^ with k < 2. In the 
easy plane case k « 1.35 from the numerical results^i 
Thus the phonons will then prempt a direct Neel- VBS 
transition and introduce a small coexistence region. It 
will thus be necessary to look at temperatures that are 
not too low to meaningfully compare with experiments. 

Barring these caveats the interesting critical phenom- 
ena discussed in this paper should be visible in a number 
of different experimental probes. Scaling forms for vari- 
ous experimentally physical observables are readily writ- 
ten down. For instance, right at the critical point the spin 
response function ^(/ci,^) near the ordering wavevector 
{Qi = (tt, tt)) will take the form 

Here ki is assumed to measure the deviation of the phys- 
ical wavevector from Qi, and T is the temperature. The 
corresponding spectral function can be directly measured 
in neutron scattering experiments. At a fixed small 
wavevector ki and temperature T, this will show sharp 
spin wave peaks as a function of frequency in the Neel 
state and similar sharp 'triplon' peaks near the spin gap 
in the VBS state. However right at the critical point, 
there will be an anamolously broad power law peak (due 
to the large rj). 
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The large value of rj will also directly manifest itself 
in NMR experiments. Indeed the nuclear spin lattice 
relaxation rate is essentially given by 

^^T [d'k hm (10.2) 

where x" is the imaginary part of the spin response func- 
tion. It is now easy to see from scaling that 1/Ti ~ 
at finite temperatures in the "quantum critical" region. 
Thus this experiment provides a direct measurement of rj. 
It is therefore an excellent way to experimentally distin- 
guish the predictions of the present paper from those of 
the earlier accepted theorjiS^ of the Neel-VBS transition 
which gives 77 « 0. 

In this context it is interesting to reconsider experi- 
ments measuring the spin- lattice relaxation rate of the 
Cu ions in the undoped and lightly doped cuprates^. 
Remarkably at high temperature the l/Ti saturates to 
a temperature and doping independent value. One sug- 
gested explanation'^'^ is that at these high temperatures 
in the undoped sample the system is in a quantum crit- 
ical regime associated with a disordering transition of 
the Neel order such that « 0. Furthermore the effects 
of doping has been suggested to only make the system 
appear closer to the critical point (at least for the high 
temperature spin physics). The results of the present pa- 
per imply that if this interpretation of the experiments in 
terms of proximity to quantum criticality is correct, then 
the corresponding transition cannot be from the Neel to 
the VBS state. 



XI. DISCUSSION 

This paper has described a variety of quantum critical 
points in two dimensions which can be understood using 
the new paradigm of 'deconfined' quantum criticality^. 
The critical point has an emergent topological conserva- 
tion law, and the critical theory is expressed most natu- 
rally in terms of fractionalized degrees of freedom. The 
order parameters characterizing the phases flanking the 
critical point emerge at large length scales as composites 
of the fractionalized modes, or their duals. These exam- 
ples clearly violate the LGW paradigm, in that the order 
parameters are not directly related to the critical modes. 

Our primary example was the Neel to VBS transi- 
tion for the = 1/2 square lattice antiferromagnet. We 
showed that a deconfined critical point scenario emerged 
in a number of tractable deformations of models appro- 
priate to describe this transition. Several results existing 
in the literaure (for instance on large-iV models) were 
shown to support this proposal when correctly inter- 
preted. We briefiy reiterate a few key physical aspects 
of this theory. First, the critical theory possesses an ex- 
tra global topoogical conservation law (associated with 
skyrmion number). It is most naturally expressed not in 
terms of the natural order parameters of either phase but 



in terms of new spin- 1/2 spinon degrees of freedom that 
are specific to the critical point. The emergence of these 
fractional spin fields at the critical fixed point manifests 
itself quantitatively in the large value of the anamoulous 
exponent 77 at the transition. The extra topological con- 
servation law is obtained because monopole events at 
which skyrmion number can change are irrelevant and 
disappear at long scales at the critical fixed point. How- 
ever they are relevant in the paramagnetic phase and lead 
to the appearnce of broken lattice symmetry. There are 
two diverging length or time scales near the critical point. 
In the paramagnetic phase the first (shorter) length is the 
spin correlation length. There is a longer length scale at 
which the VBS order gets pinned. 

We also considered a number of other critical points in 
this paper. The transition between VBS and spin liquid 
states was discussed in Section I Villi and described by 
critical theory closely related to that for the SJ mod- 
els. The superfiuid-insulator transition of bosons at 
half-filling on the square lattice was considered in Sec- 
tion llX AI the insulator exhibits bond-density- wave order 
and the theory for the critical point (Eq. had been 

obtained earlier^. Here we provided a physical reinter- 
pretation of this theory, and showed how it could also be 
understood as a deconfined QCP. 

Overall our work shows that such deconfined quantum 
criticality is quite common in two dimensional systems 
with a spin-1/2 moment per unit cell. This leads us to 
suspect that the scope for finding deconfined QCP's in 
other correlated electron systems, including those with 
fermionic excitations, is bright. These QCP's naturally 
have large anomalous dimensions for observable order 
parameters^, and so hold the prospect of explaining a 
variety of experimental puzzles. 
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APPENDIX A: BERRY PHASES IN THE SJ 
MODEL 

The non-linear sigma model representation in Eq. I|3.1|) 
associates the Berry phases with a summation over the 
individual Berry phases of each spin. Each such contribu- 
tion measures the area enclosed by the path of the spin on 
the unit sphere, and this is represented by Sb in Eq. (|3.1|l . 
Upon transforming to the z variables via Eq. H1.5|l , there 
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is a simple way of measuring this axeaSSiS^i^: it is the 
Polyakov loop of the U(l) gauge field of the CP^ model. 
This connection suggests that the Sb in Eq. (j3.13|l should 
be replaced by 



S'b = 



(Al) 



where is the cubic lattice, r-independent representa- 
tion of the square lattice sublattice staggering factor er- 
One can now consider a 'modified' SJ model'^^ with ac- 
tion 5^ J ^ + Sa + S'g defined by Eqs. |TTT|l . (|TT^ 
and ljAl|l . This appendix will argue that the properties 
of Ss] and S'gj are very similar, and universal features 
are identical. 

First, we show that for i = (this is deep in the VBS 
phase), the two theories are, in fact, exactly equivalent. 
In this limit, we can proceed with a duality mapping as 
in Section [iVAl and obtain a dual representation of Ssj 
which is Eq. (|4.9|) with A = 0: 



5= 2^(A(x-K.?))' 



(A2) 



Proceeding with the analogous duality transformation to 



5gj, we find instead 



(A3) 



where Bq is a fixed integer-valued field on the links of 
the dual lattice chosen so that 



A X Bo = 



(A4) 



where r is a unit vector in the r direction. A convenient 
choice for Bq is shown in FigEt.. Now note that we can 
write 



Bo = A?9 + A X /3 



(A5) 

with -dr defined below Eq. (|4.4|l and shown in Fig^ /3 is a 
fixed vector field on the links of the dual lattice with only 
its temporal components non-zero as shown in Fig ^p. 
It is now evident that Eqs. I|A3|I and l|A5|l are exactly 
equivalent to Eq. (|A2|) . as the couplings between (3 and 
X, "& vanish. 

Moving to the general case with < 7^ 0, let us exam- 
ine the fate of the modified SJ model for TV = 1 as in 
Section ITVAI In this case, Eq. (|4.9|l is replaced by 



S' 



—.(Ax A? 

2t ^ ' 



1 

2K 



(A(x + i9) + A + A X /3)^ 



(A6) 

It is now not difhcult to show that the additional term 
associated with /3 above makes little difference to the uni- 
versal properties of theory: integrating out the massive 
A modes is now a little more involved, but the final the- 
ory for X has the same structure as that in Section llVAl 
Further details may be found in Ref. I20I 

Similar comments apply to the modified SJ model at 
N = 2, with easy plane anisotropy, as discussed in Ap- 
pendix 



X X X ^ X 

X i"^^ X 
• • • • • 

(a) 

+1/8 -1/8 +1/8 -1/8 

XXX 
-1/8 +1/8 -1/8 +1/8 

XXX 
+1/8 -1/8 +1/8 -1/8 
(b) 

FIG. 6: Specification of the non-zero values of the fixed fields 
(a) Bo and (b) f3. The circles are the sites of the direct lattice, 
j, while the crosses are the sites of the dual lattice, J; the latter 
are also ofi'set by half a lattice spacing in the direction out of 
the paper (the /i = r direction). The Boi_i are all zero for 
fi = T,x, while the only non-zero values of Boy are shown in 
(a). Only the fi = t components of /3p are non-zero, and these 
are shown in (b). 



APPENDIX B: DUALITY TRANSFORMATION 
WITH EASY PLANE ANISOTROPY 

The duality transformation for the SJ model at = 2 
in the easy plane limit proceeds very similarly to that 
of the A^ = 1 case discussed in Section II V Al We begin 
by rewriting the 'boson hopping' term in Eq. H5.5f) in a 
Villain approximation: 



E 



-ijc • (A(/)„ - a) 



(Bl) 



where a = 1,2 labels the two species of bosons and jo, 
are the corresponding integer valued currents defined on 
the links of the square lattice. Proceeding as before in 
Section IIV Al integrating out the (pa fields leads to the 
current conservation conditions 



A-j„ = 0. 

which can be solved by writing 

= A X Aa 



(B2) 



(B3) 



with Aa integer fields. We treat the gauge field kinetic 
energy term as in Section flV Al bv first decoupling it with 
the b field, and then summing over the integer q to obtain 



b A?9 = B 



(B4) 
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with B an integer field. Integrating over the gauge field 
a now replaces Eq. (|4.7II by 



A X B = ji +j2 
This may be solved by writing 

B = Ai + A2 + Ax 
with X ^-ii integer. The action then reads 



(B5) 



(B6) 



S 



E 



i(A X Ai)^ + i(A X A2)^ 
+ ^(A(x + i?) + Ai+A2f 



(B7) 



As in Eq. H4.10|l we may soften the integer constraints 
on Aq, , X by adding the terms 

—t cos(27rAi ) — t cos(27r A2) — ^ A„ cos(27rnx) (B8) 



with n running over all positive integers. Now we can 
shift X X — X + ^- Then we can put 27rx = O1 — 62, and 
integrate over both phase fields, Oa, leaving the partition 
function unchanged up to an overall multiplicative con- 
stant. Upon shifting the two fields Ai Ai — A.9i/2tt 
and A2 ^ A2 + A02/27r, the last term in Eq. (jB7|l takes 
the form (Ai + A2)^. We can then define A+ = Ai + A2 
and A = 7r(Ai — A2), and integrate out the massive field 
A+ . Up to irrelevant terms we thereby obtain for the full 
action: 



— (Ax A)^ - icos(A6'i - A) 



tcos(A6l2 - A) 



A„ cos(n(6li -02- 2n'&)) 



(B9) 



Once again as e oscillates on four sublattices, for 
smooth variations of O12, the lowest value of n that sur- 
vives is n = 4. We may therefore replace S\ by 

Sx = -J2[Xcos{A{0i-02))] (BIO) 

with A = A4. The resulting action is then a 'hard-spin' 
version of the action Sz in Eq. lESll of Section El with 
the identification of the vortex operators, ^pa ~ e'^° . 



APPENDIX C: ESTIMATE OF MONOPOLE 
SCALING DIMENSION 

We can ask about the answer for all these scaling di- 
mensions that would be obtained by first integrating out 
the matter fields, truncating the resulting gauge action 



to quadratic order, and using that theory to calculate the 
scaling dimension. This would be an estimate, albeit un- 
controlled; one might hope that it will correctly capture 
the trends. 

Anyway, it turns out that this can be done without 
serious calculation. First, note that in this procedure the 
answer depends only on the universal conductivity of the 
matter fields at the transition ignoring all coupling to 
the gauge field. The higher this universal conductivity 
the higher the instanton anomalous dimension. 

The simplest case is the iV = 1 SJ model. Here the 
relevant universal conductivity is that of a single boson 
species. We know that at this transition the four in- 
stanton operator is irrelevant. Now in the N — 2 cases 
with either easy plane or full SU(2) symmetry, it is clear 
that the universal conductivity will only be larger than 
at iV = 1 . Thus we should expect a higher anamolous di- 
mension. This would then predict irrelevance of the four 
instanton term with or without easy plane anisotropy for 
the physical case of = 2 in agreement with other ex- 
pectations. 



APPENDIX D: SJ MODELS IN ONE DIMENSION 

There is a close and useful analogy between some of 
the phenomena explored in this paper and correspond- 
ing ones in one spatial dimension. Specifically consider a 
one dimensional spin- 1/2 magnet in the presence of some 
easy plane anisotropy. The analog of the Neel phase in 
rf = 1 is a phase with power law correlators for the stag- 
gered XY magnetization. This phase may be described 
as a Luttinger liquid. There is a direct second order tran- 
sition between this phase and a VBS phase where there 
is spontaneous dimerization of the spin chain. 

A useful theoretical description of this transition is ob- 
tained by focusing on vortices in the space-time configu- 
ration of the staggered XY order parameter field. From 
a quantum mechanical point of view such vortices cor- 
respond to phase-slip (or instanton) events. It is well- 
known that a 27r phase slip event carries a momentum 
TT, and hence is not allowed as a term in the Hamilto- 
nian. (In an equivalent description of this Luttinger liq- 
uid phase in terms of interacting spinless fermions, these 
2tt phase slips correspond to interchange of left and right 
movers). In the VBS phase these phase slip events have 
proliferated. Indeed it is precisely the tt momentum that 
is carried by the 27r phase slip that is responsible for the 
broken translation symmetry of the VBS. A convenient 
order parameter for the VBS phase is therefore provided 
by the 2tt phase slip operator. 

Though 2tt phase slip terms are not allowed in the 
Hamiltonian 47r phase slips (which carry zero crystal mo- 
mentum) are clearly allowed. The transition from the 
Luttinger liquid to the VBS is driven by the proliferation 
of these An phase slips. These 'doubled' phase slips are 
irrelevant throughout the Luttinger liquid phase and are 
marginally irrelevant at the critical fixed point. 
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The analogy with the two dimensional situations con- 
sidered in this paper is now clear. The 27r phase slip is 
the direct analog of the skyrmion tunneling {i.e single in- 
stanton) event. In both d — 1,2 the VBS phase is under- 
stood as a condensate of the appropriate single instanton 
event. In d = 1 the transition is driven by doubled in- 
stanton events (similar to the quadrupling of instantons 
in two dimensions) which stay irrelevant at the critical 
fixed point. At a formal level it is possible to construct 
an appropriate 'SJ' model that correctly describes the 
transition even in d = 1. 

A complete presentation closely related to the discus- 
sion in this appendix appears in Ref. i38L 



APPENDIX E: DIRECT DERIVATION OF DUAL 
MERON ACTION 

We pass from Eq. H9.1|l to the analog of a CP^ repre- 
sentation by letting 



'(<?>lr— <?>2r 

nir — n2r + 1 



nir + n2r 



(El) 
(E2) 
(E3) 



Here 6i_2 — e"^^'^ represent charge-±l/2 bosonic opera- 
tors and ni^2 are the corresponding boson numbers, and 
Er was defined in Eq. (|1.3|l . Note that bi^2 are not canon- 
ical Bose operators, and the relevant commutation rela- 
tions here are [ni, (pi] = —i and [n2, ^2] = —i- The eigen- 
values of ni^2 are integers which run from —00 to 00. As 
is usual, there is a gauge redundancy associated with an 
arbitrary choice of the local phase of the 61^2 fields. The 
last equation is a constraint that requires the total num- 
ber of both species of bosons to be fixed at -1-1 on the 
A sublattice and —1 on the B sublattice. The left-hand 
side of this constraint equation is precisely the generator 
of the local gauge transformation. We have chosen to 
stagger this gauge charge on the two sublattices. 

The Hamiltonian in Eq. (|9.1|l is readily rewritten in 
terms of these new variables: 



H 
Hu 



Hu + Ht 
U 



(E4) 



"2 X! ~ <^rnof + {n2r - er"-o)^ (E5) 



(rr') 



(E6) 



We have introduced a term proportional to uq which de- 
scribes a 'chemical potential' for the total on-site gauge 
charge. As the total gauge charge is fixed to ±1 on each 
site, this addition is completely innocuous (for any value 
of no). Later we will choose hq appropriately to ensure 
that the ni_2 fields have zero mean value. While this step 
is not necessary, it is convenient, and will be commented 
upon further at the appropriate point. 



Now we proceed to a path integral representation to 
write 

S = Su + Sr + Sao + 

5(7 = ^ fdr ^ ^ ^{nir - e^no)^ -|- {n2r - CrUof 

Sr — ^ / dr [inirdrCplr + in2rdr(f>2r] 
r '' 

Sao = X! dr [iao {nir + n2r - er)] 

r •' 

St - [ drHt (E7) 



As usual, it is assumed that there is a fine discretiza- 
tion of the imaginary time index t, and that the integer- 
valued boson numbers nir, «2r five on the temporal links 
at each spatial point. The 'gauge' constraint is imposed 
by means of a Lagrange multiplier field ao which will be 
interpreted as the time-component of a gauge field. We 
now proceed as is usual in slave particle theories of corre- 
lated systems. We decouple the interactions in Ht using 
a complex auxiliary field Xrr' defined on each spatial link 
to write 



.-St 



n / dXr 

(rr') 

2 

+Xrr' ^ I 



exp 



\Xr 



(E8) 



with Xrr' = Xr'r- The fluctuations in the amplitude of 
the X field are expected to be innocuous. Hence we will 
write 



Xr 



Xoe 



(E9) 



with xq a constant that simply renormalizes the boson 
hopping amplitude. As usual the will be in- 

terpreted as the spatial component of a gauge field. The 
full action now is invariant under the gauge transforma- 
tion 



ao 



g«7r(T-)gi0„r(-r) 



dr 

-Irir) +'-fr'{T), 



(ElO) 
(Ell) 
(E12) 



where a = 1,2. To examine universal critical properties 
near the phase transitions of interest it is legitimate to 
add various possible local terms that are consistent with 
the global symmetries and gauge structure of the action. 
It is particularly useful to add a 'kinetic energy' term for 
the gauge fields on all plaquettes (spatial and space-time) 
in Villain form: 



Ce = -Ef + lE, {dra, ~ A,ao) 



-B^ + iB {e^jAiaJ 



(E13) 
(E14) 
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Here Ei is an integer-valued 'electric field' defined on 
the spatial links at each time-slice extend over the 
spatial co-ordinates x, y) and B is the corresponding 
integer-valued magnetic field on a spatial plaquette. We 
have introduced the vector notation Ui — {ax,ay) — 
(ap^P^x 1 o.r,f+y) ■ For simplicity we have chosen the same 
constant u multiplying the Ef and terms. The origi- 
nal microscopic action is formally obtained in the large-u 
limit {u — > oo). 

To proceed it is first useful to note that the background 
gauge charge present in this formulation will lead to a 
background electric field about which the true electric 
field will actually fluctuate. It will be convenient to in- 
corporate this effect by finding a suitable mean field for 
the various fields in the action. Consider a mean field 
(saddle point) of the action where the non-zero expecta- 
tion values are 

(m) = (n2) = n 
(ao) = ao 

{E,) = Em- (E15) 

The saddle point equations are obtained by varying the 
action with respect to these fields: 

uEio = iAitto 
AiEiO = -2n + er 
Uin — ernQ) = — (E16) 

We now use our freedom in choosing the constant no to 
set 



so that n — Q. One may worry that this special choice of 
no might indicate some non-generic nature of the result- 
ing theory. We note, however, that qualitatively identical 
results are obtained for any other choice of no - in the 
dual action with such a choice the merons see some non- 
zero but spatially oscillating flux. Having a zero spatial 
average, this flux does not qualitatively effect the low en- 
ergy (extended) meron states. The above choice simply 
renders the low energy behavior more transparent. We 
then have 

A,S,o = (E18) 

Note also that i?io is the gradient of a potential deter- 
mined by ao. These conditions determine Ei^ (the back- 
ground electric field) uniquely to have the value 1/4 ori- 
ented from the AtoB sublattice. 

We may now examine the full theory by first shifting 
Oo = Oo + . Straightforward manipulation shows that 
the electric field now fluctuates about a background value 
EiQ so that the Ei dependent terms in the action read 

= ^ {E, - E,^f -I- lE, {dra, - A,ao) (E19) 

All other terms remain unchanged (after the replacing oq 
with (Sao). We may now dualize this action as in previ- 
ous sections to directly derive the dual meron action of 
Section El 



APPENDIX F: BREAKDOWN OF THE 
'SCREENING ARGUMENT' IN THE 
MONOPOLE GAS 

In this Appendix we will consider a simple toy model 
of a compact U(l) gauge theory without Berry phases 
which can be shown to possess a deconfined critical point. 
This will enable us to understand clearly the claim of 
Section IIVBI that the specific monopole gases that ob- 
tain at the critical points studied in this paper evade the 
general monopole screening argumenta^Siffi for a three 
dimensional Coulomb gas with logarithmic interactions. 

We consider a model of charge n bosons (n > 4) cou- 
pled to a compact U(l) gauge fleld in Z? = 2 + 1 dimen- 
sions with Euclidean action 

5 = -J^cos(A(/)- na) - fi:^cos(A X a). (Fl) 

p 

Here e*"^' represents a boson field on the sites i of a three 
dimensional cubic lattice. The sum in the first term is 
over the links of the lattice, while that in the second 
term is over the plaquettes P. The field a^- = a^ -I- 2tt 
is a compact U(l) gauge field, and the integer n is the 
charge of the boson. This model has two phases. For 
large J, there is a Higgs phase where the boson field 
has 'condensed'. Following the arguments of Ref . IbtI the 
effective theory of this phase is a Z„ gauge theory in its 
deconfined phase in 2 4- 1 dimensions. The excitations 
in this phase are stable vortices that carry flux 27rq/n, 

for q = 1, n — 1, of the gauge field a. For small J, 

on the other hand, there is a different phase which is 
associated with confinement of the U(l) gauge theory. In 
particular, the Zn vortices that appear in the Higgs phase 
are no long er present in the spectrum. As also argued 
in Rcf. 1671 the transition between these two phases is 
actually described by that in a Zn gauge theory. This 
latter theory is dual to the global Zn clock model - for n > 
4 the clock anisotropy is irrelevant and the universality 
class is 3D XY. 

The n = 1 case of Eq. ljFl|l was considered in Ref. |2^ 
For this case, the 'clock' anisotropy is strongly relevant 
(it rounds out the transition into a crossover), and the 
physics is very different from n > 4. 

Formally, the action above is readily dualized (as in 
many of the other examples discussed at length in previ- 
ous sections). The dual action takes the form 

5duai = 51 [cos(2^Ax) - A cos(27rnx)] (F2) 

and has the expected structure of a global XY model with 
n-fold anisotropy. 

It is useful, for our purposes, to have a physical in- 
terpretation of these results. In the Higgs phase, the 
vorticity of the (j) field is quantized in units of 27r/n (as is 
natural for a charge n condensate). The presence of in- 
stantons implies that the total flux can change by 2tt, so 
that n of these vortices can appear or disappear together. 
Thus the vortices only carry a Zn quantum number. The 
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dual description focuses on these vortices. Without 
instantons, the I-k jn vortex is globally conserved and its 
physics is described by a global XY model (this is just 
the duality in Ref . Issf) . The presence of instantons leads 
to the n-fold anisotropy (the A term in Eq. ljF2p ) for 
this global XY model, leading to the global Z„ model. 
Thus the irrelevance of the n-fold anisotropy, for n > 4, 
should be interpreted as the irrelevance of instantons at 
the transition. Indeed the XY universality class, is the 
exact dual of the condensation transition of the charge n 
boson coupled to a non-compact U(l) gauge fieldiS 

Now let us analyze the transition in the 'RPA' ap- 
proach outlined in Section llVBI The transition is associ- 
ated with the condensation of the e*"^ field. We therefore 
integrate out this field in the presence of a non-trivial 
gauge potential, and truncate the resulting gauge action 
to quadratic order (initially ignoring instantons). The 
result is, as in Eq. (|4.12l) . 

^aon'\K\\M^)\' + - (F3) 

Here ctq is a universal constant, and slt is the trans- 
verse component of the gauge field. The Maxwell term 
present in the bare action is less important at long dis- 
tances than the term in the action displayed above, and 
we have dropped it. We now examine the stability of Sc 
to instantons. First, we dualize Sq to obtain22iSS 

— \X{K)\^ - Y,\cos{2T:nx). (F4) 

The last term represents instanton events. Note the 
in the first term. Now the A term will in general gen- 
erate a term in the Gaussian x action which will 
then eventually make instantons relevant. A fine-tuning 
is required22i2£ to prevent the generation of the term 



in this argument, and for the present model we can now 
easily see that this fine tuning is automatic at the critical 
point of the gauge theory in Eq. (|F1(I . 

The key is to note that the logarithmic interaction be- 
tween the monopoles is equivalent to the statement that 
the correlators of e^*'^"-^ decay as a power law (at the 
fixed point without monopoles). In the theory of the 
monopoles in Eq. ljF2|l . the e^^^^ field is at the critical 
point of the 3D XY model. The RPA theory approx- 
imates this non-trivial interacting critical theory by an 
equivalent Gaussian theory, which also happens to give 
power law correlations for the e^™^ field (this is a prop- 
erty of the iiT^lxp form). In the full theory in Eq. I|F2|I . 
the screening of monopole interactions is associated with 
corrections higher order in A. In the context of conformal 
perturbation theory about the critical point of Eq. I|F2|) . 
however, it is clear that these higher order effects in A 
actually represent shifts in the position of the critical 
point, and not any changes in the scaling dimensions of 
operators. 

Hence for the Neel-VBS transition we conclude that 
the 'naive' computation of monopole scaling dimensions 
in the large N limiliS is actually correct, and that we 
should neglect screening between multiple monopoles in 
determining this scaling dimension. The latter effects are 
more correctly accounted for by shifting the position of 
the critical point. 

Note that for the n — \ case of Eq. (jFl|) considered 
in Ref. 28, computation of the scaling dimension of the 
monopole operator using Eq. I)F2|I shows that monopoles 
are relevant at the critical point. Indeed, they round out 
the transition to a crossover, and the monopoles are al- 
ways in a screened plasma phase. So the conclusions of 
Refs. I29ll3(1 for this case are correct, but not for com- 
pletely sound reasons. 
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